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PREFACE 


The  Seventeenth  DoD  Tri-Service  Review 
Conference  on  Atmospheric  Transmission  Models  was  held 
at  the  Geophysics  Directorate  of  the  USAF  Phillips 
Laboratory  at  Hanscom  AFB,  Massachusetts  on  7-8  June 
1994. 


More  than  190  authors  and  participants 
presented  45  papers  in  the  field  of  the  coupling  and 
interaction  of  optical  and  atmospheric  phenomena. 
Sessions  were  held  on  Atmospheric  Radiation  Codes; 
Aerosols, Clouds  and  Scattering;  Structure  and 
Turbulence;  High-Altitude  Models  and  Measurements; 
Environmental  Applications:  a  Poster  Session  featuring 
other  topics  was  also  held. 

These  volumes  include  the  abstracts  and  hard¬ 
copy  and  any  other  material  provided  by  the  authors 
for  their  presentations.  An  Author  Index  and  a  partly 
overlapping  attendee's  list  are  also  provided. 

The  abstracts  are  grouped  at  the  back  of  the 
second  volume  on  pages  A1-A45  for  ease  of  reference 
and  perusal.  The  author  list  includes  both  abstract 
(A#)  and  paper  references.  Because  of  the  delay 
between  the  submission  of  the  abstract  and  the  hard¬ 
copy  there  is  not  always  an  exact  correspondence 
between  the  two. 


jsail  P.  Anderson 
Simulation  Branch 
Optical  Environment 
Division 
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17th  ANNUAL  TMOSPHERIC  TRANSMISSION  MODELING 
ANNUAL  REVIEW  CONFERENCE 

7  JUNE  1994 

GEOPHYSICS  DIRECTORATE/PHILLIPS  LABORATORY 


FAS  CODE:  Validation  and  Applications 


G.P.  Anderson,  J.H.  Chetwynd,  J.  Wang* 
Geophysics  Directorate/PL 
*NRC  Fellow  (on  leave  from  Univ.  of  Michigan) 


J.-M.  Theriault 

DREV/Defence  Research  Establishment  Valcartier 


L.W.  Abreu 
ONTAR,  Inc. 

S.A.  Clough,  J.-L.  Moncet 
Atmospheric  and  Environmental  Research,  Inc. 
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FASCOD3:  Its  Future 


Who  are  we? 

Gail  Anderson,  Jim  Chetwynd,  Steve  Miller,  Jinxue  Wang 
PL/GPOS 

What  are  we  doing?  (FY94/95) 

VALIDATION 

GEOMETRY 

LASER/LIDAR  Specificity 
Inversion  Algorithm  Development 
Collaboration  with  DOE’s  LBLRTM 

How  are  we  doing  it? 

LBLRTM  =  Line-by-Line  Radiative  Transfer  Model 

Predicated  on  FASCODE;  Funded  for  last  4  yrs  by  DOE 
Responsible  AER  Authors:  Clough,  Worsham,  Moncet 
EASE:  FAScode  for  the  Environment;  Moncet,  Gallery  (AER) 
and  Wang  (GPOS/NRC) 

Combination  of  best  in  both  LBL  Codes: 

non-LTE,  geometry,  &  lidar  applications 
from  FASCODE 

vectorization,  radiance  algortithm,  optimization 
from  LBLRTM 

Inversion:  Miller  (GPOS)  and  Moncet  (AER); 

also  Theriault  (DREV)  and  Wise. 

Nadir:  Univ.  of  Wise.  HIS  Bomem  ER-2  Spectrometer 
Up-Looking:  Wise.  AERI  &  HIS;  Miller 
Limb:  Nadile  (CIRRIS-IA);  O’Neil  (MSX) 

Validation: 

IR:  WMO  ICRCCM  (Ellingson)  IR  by  Wang  (GPOS) 
mm:  WMO  ICRCCM  (Westwater)  mm  by  Chetwynd,  Hoke 
(GPOS)  &  Clough  (AER) 
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FASCODE:  Fast  Atmospheric  Signature  Code 
Version  3  -  ^ 


AFGL:  F.X.  Kneizys 

S.A.  Clough^ 

G.P.Anderson 
J.H.  Chetwynd 
L.W.  Abreu 
W.O.  Gallery^ 

L. S.  Rothman 

M. L.  Hoke 
E.P.  Shettle^ 
R.D.  Worsham^ 


1  -  Primary 

I 

_j 

-  UV,  Constituent  Profiles 

-  Programming,  Validation 

-  LOWTRAN?  Compatibility 

-  Geometry 

-  Line  Atlas 

-  Line  Coupling 

-  Aerosols/Hydrometeors 

-  X-sections,  etc. 


1  currently  at  Atmospheric  and  Environmental  Research,  Inc. 
^  currently  at  Naval  Research  Laboratory 
^  at  Atmospheric  and  Environmental  Research,  Inc. 
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FASCODE:  Fast  Atmospheric  Signature  Code 
Version  2  and  Version  3-/9 


Visidyne: 

H  J.P.  Smith 

D  J.  Dube 

M.E.  Gardner 

-  contributions  to  FASCODIC  (1978) 

T.C.  Degges 

-  NLTE  theory  (1977, 1985) 

Sonicraft: 

W.L.  Ridgway 

R.A.  Moose 

A.C  Cogley 

-  contributions  to  FASCOD2  (1985) 

AER,  Inc.: 

R.G.  Isaacs 

-  multiple  scattering  (1987) 

R.D.  Worsham 

-  programming  contributions  to 

FASCOD3  &  MS  (1988-91) 

S.A.  Clough 

-  radiance  algorithm,  etc.  (1990-1991) 
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HISTORY 


DoD  Plan  for  Atmospheric 
Transmission  Research  and  Development 

AIR  FORCE 

o  Maintain  DoD  Standard  Atmospheric  Optical/IR  models: 

(LOWTRAN),  MODTRAN,  FASCODE,  HITRAN  Database 
o  Publish  and  Brief  Model  Updates 
o  Conduct  Annual  Tri-Service  Review 
o  Measure  and  Model  Propagation  Effects  of  the 
Free  Atmsophere 

ARMY 

o  Study  Battlefield  Conditions 

o  Develop  Models  of  Dust,  Smoke,  Chemicals,  Propagation, 
and  Diffusion  Effects 

NAVY 

o  Develop  Models  for  Marine  Environment 
o  Measure/Model  Atmospheric  Propagation 
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The  Problem:  RADIATIVE  TRANSFER 


1.  The  Atmosphere  as  a  contaminant  for  E/0  Systems 

2.  The  Atmosphere  as  a  signature  source  for  natural  variability 

Solutions:  DEFINITIONS 

1 .  State  variables  (T,  p,  Cld,  Aer)  along  line-of-sight 

2.  Spectral  Characteristics  of  the  Path  Variables 

3.  Viewing  Geometry 

4.  E/0  System  Characteristics  (Spectral  Range  &  Resolution, 

Platform,  Objective) 

Solutions:  OPTIMIZATION 

1 .  Efficient  Mathematical  Algorithms  (Line-by-Line) 

2.  Accurate  Band  Model  Options 

3.  User  Friendly 

4.  Validation/Documentation 

Solutions:  DATA  ANALYSIS 

1 .  Information  Theory 

2.  Inversion  Algorithm  Development 

3.  Ground  Truth 

4.  Validation  and  Error  Estimation 
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DEFINITIONS 


Ki  =  absorption  cross  section,  related  to  molecular  properties, 
pressure  (p),  temperature  (6) 

rji  =  column  amount  of  absorbing  (i’th)  species  =  J  n^ds 
ds  =  path  increment;  n„  =  volumn  density 
Ti  =  optical  depth  =  ic;  17; 

Tj  =  transmittance  =  exp  (-T;) 

Tmoi  “  molecular  transmittance  =  t  T;  =  T,  T2  T3  T4..... 
Tt  =  total  transmittance  =  Taerosoi 

B(0)  =  Planck  Function  for  temperature  0 
=  Solar  (Lunar)  Source  Function 
^  =  Non-LTE  Source  Term 
W  =  Weighting  Function  =  (cTT-r/ds) 

OPERATIVE  EQUATION: 

=  Thermal  Radiance  =  {  B(0)  dTj  =  {  B(0)  (dTj/ds)  ds 
AND  FOR  A  SINGLE  LAYER: 

^  =  {  B(I9)  dTj  =  B(0)  [1-Tt] 
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PAGE  2 


and,  finally,  combining  thermal,  solar,  and  non-LTE  sources 
with  multiple  scattering, 

one  can  replace  the  Planck  source  function  with  a  more  general 
source  function: 

B(^)  /(r,r) 

7(r,f)  =  general  source  function  dependent  on  optical  depth  (t)  and 
viewing  geometry  (f  =  zenith  &  azimuth  cosines) 

=  T^T(f)  p  [SOLAR] 

+  [l-ojJ  Bid)  [THERMAL] 

+  ^  [N-LTE] 

+  Jj^s  [MUL.SCAT.] 


where: 

=  single  scattering  albedo  =  /(r^+rj 

T(f)  =  transmittance  from  top  of  atmosphere  to  layer 
p  =  scattering  phase  function 

Jms  —  for  both  solar  and  thermal  terms  (complicated) 
and  other  terms  are  as  previously  defined!! 


O 
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*  VALIDATION  for  Signatures 

1.  targets /pollutants  against  backgrounds- 

estimate  of  variability 

estimate  of  go /no go 

estimate  of  true/false;  reduce  false  ID 

incorporate  "target  signatures " 

2.  atmospheric  specification 

local  and/or  global  "weather" 
atmospheric  contaminants  become  signatures 

3.  higher  resolution  spectroscopy  with  high  spectral 

accuracy  may  minimize  confusion,  increasing 
signal/noise  for  both  contaminants  and  state 
specification. 

*  REALTIME  ANALYSIS  with  speed 

1.  Line-by-line  (LBL)  codes  and  their  derivatives  work  - 

"exact" physics  must  remain  state-of-the-art; 

Because  physicists  drive  technology  to  higher  resolution: 
lidar/laser  applications 
intetferometers  at  <.01cm-l  resolution 
remote  sensing  issues:  line  coupling,  mm-wave 
new/old/unsettled  issues: 

CFC  x-sections,  H20  continua,  C02  x  factor 
NLTE  issues 

However: 

"LBL"  is  slow,  slow,  slow!!! 

2.  Band  Models  (pragmatic , expedient  parameterizations)  trail 
"LBL"  in  "state-of-the-art", 

but  they  are  the  codes  of  the  NOW  and  FUTURE  for 
issues  and  analyses  with  compatible  spectral  resolution! 
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Figure  1 
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WAVENUMBER(cm^) 


SPECTRE  EXPERIMENT:  case2 


(,.ujo-Jts-,ui/MUJ)  30NViaVU 


Figure  2 
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WAVENUMBER(cmO 


Figure  3 
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SPECTRE  EXPERIMENT;  case4 


(,.mo-jjs-,LU/Miu)  aONViaVU 


Figure  4 
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WAVENUMBER(cm’) 


BEUGHTNESS  TEMPERATURE  (K) 


FASCOD3  CALCULATED  BT  CO 


Fig.  5.  Correlation  plots  of  FASCOD3  calculated  brightness  temperatures  and 
measured  brightness  temperatures  supplied  by  E.  Westwaier. 
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FASCODE  Introduction 


FASCODE  -  FASCOD3P  (lQ-1992) 
Fast  Atmospheric  Signature  CODE 
Version  3  -  Prelimary 


POINTS  OF  CONTACT: 

Gail  P.  Anderson/James  H.  Chetwynd 
Simulation  Branch 
Optical  Environment  Division 
Geophysics  Directorate/PL 
29  Randolph  Road 
Hanscom  AFB,  MA  01731-3010 
617-377-2335/2613;  FAX  xxx-8900. 
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FASCODE 


OBJECTIVE:  To  provide: 

o  "exact"  spectroscopic  calculations  of: 

molecular  radiative  transfer  properties 
Voigt  line  shape 

LTE  and  non-local  thermodynamic  equilibrium  (non-LTE) 
line-by-line  (LBL)  code 

heavy  molecule  x-sections  (i.e.  chloro-fluorocarbons,  etc.) 
laboratory  and  atmospheric  simulations 
provisions  for: 

Non-molecular  (clouds,  aerosols,  geometry,  etc.) 
from  LOWTRAN/MODTRAN 
altitude  range:  open 
lines-of-sight:  open 

spectroscopic  constants  from  HITRAN92  and/or  other  sources 
x-sections  from  HITRAN92  and/or  other  sources 
non-LTE  populations  from  SHARC 
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FASCODE 


TECHNICAL  DESCRIPTION: 

FASCODE  calculates  atmospheric  radiance  (3?)  and  transmittance  (T) 
for  the  following  conditions: 

full  Voigt  spectral  resolution  from  0  to  50000  cm-1 
"exact"  line-by-line  radiative  transfer  in  UV/Vis/IR 
optimized  line  shape  and  layering 
scattering:  Rayleigh,  Mie,  single,  multiple 
6  climatological  descriptions:  tropical,  mid-lat  sum/win, 
subarctic  sum/win,  and  US  Std.  for 
6  atmospheric  gases:  H20,  C02,  03,  N20,  CO,  CH4,  plus 
a  "default"  profile  for  each  HITRAN  gas  (up  to  35  species) 

10  "greenhouse"  (CFC)  gas  cross-sections  and  profiles 
aerosol  profiles: 

tropospheric:  rural,  urban,  desert.  Navy,  fogs, 
stratospheric:  background,  volcanic  (background, 
aged,  high,  fresh,  and  extreme) 
clouds/rain:  cumulus,  altostratus,  stratus,  strato  cumulus, 
nimbostratus,  cirrus  (standard,  subvisual  and  NOAA) 
geometric  lines-of-sight:  HI  (observer  location)  to  H2 
(end  of  path)  with  HI  or  H2  =  surface,  space,  or 
anyplace  within. 

radiance  sources:  LTE,  non-LTE,  thermal  and  surface  radiation 

primitive  target  specification 

laser  applications 

ground  emittance  and  reflectance 

user  profile  specifications 

instrument  specifications:  plot,  scan,  filter  options 
weighting  functions,  primitive  flux  calculations 
external  data  sources:  HITRAN  and  SHARC 


19 


FASCODE 


BACKGROUND: 
o  Requirement: 

an  exact  spectral  simulation  code,  containing  the  true 
physical  equations  of  radiative  transfer. 

DoD  standard  for  judging  the  accuracy  of 
more  pragmatic  and/or  efficient  algorithms 

**  need  cannot  be  overstated  ** 
provides  a  window  upon  the  actual  workings 
of  the  atmosphere  (and  lab) 


o  Penalty: 

FASCODE  is  ponderous  and  time-consuming 
no  solar  capability 

large  storage  requirements  for  some  paths 

0  Code  development/maintenance  is  mandated  by: 

DoD  1978/83  Atmospheric  Transmission  Plan 
no  single  "sponsor"  for  the  code. 


FASCODE 


ASSUMPTIONS/CAVEATS: 
o  Accuracy: 

FASCODE  is  spectroscopicly  a  1-3%  code,  with  caveats 

line  shape  description  is  generally  within  1  % 

end-to-end  accuracy: 

spectroscopic  line  parameters  (1-3%) 
scattering  algorithms  (10%) 

LTE  field  descriptors,  p,  T,  N(i),(all  at  5-10%) 
non-LTE  field  descriptors,  excited  state  populations 
from  SHARC  (approximately  10%) 
governing  physical  equations,  molecular  (1-5%), 
particulate  (5-20%) 


NOTE:  Discrepancies  with  measurements  can  often  be  used  to  infer 
the  state  of  the  atmosphere  through  inversion  or  to  define  the  need 
for  improved  spectroscopic  identifications. 

In  summary:  altitude  limits  0-4- 120km  (LTE,  non-LTE) 

spectral  limits  0-50000cm-l  (A  500cm- 1) 

accuracy  1-3%  for  spectroscopy 

field  descriptors 

LTE&NLTE  5-15% 

applications  thermal  radiance 

transmittance 
contrast  studies 
surveillance 
plumes 
laser/lidar 

high  spectral  resolution 
microwave 


FASCODE 


CURRENT  STATUS: 

o  FASCOD3P  access: 

-  available  in  FORTRAN  for  mainframe  or  workstation 

-  implemented  on  Cyber,  Cray,  Apollo,  VAX,  and  SUN 

-  NOS,  VMS,  and  UNIXAJLTRIX  operating  systems 

-  transfer  medium:  9-track  magnetic  tape 

successful  in-house  PC  version  is  under  development 


REMINDERS: 
o  negatives: 

-  not  particularly  user-friendly 

-  computer  time  and  storage  intensive 

-  no  solar  capabilty 

o  balance  between  "first  principles"  scientific  algorithms 

versus  more  expedient  and  rapid  codes  (such  as  MODTRAN,  LOWTRAN, 
SHARC,  and  their  derivatives) 

o  singular  capabilities: 

-  15  micron  C02  line  coupling 

-  02  mm  line  coupling 

-  chloro-fluorocarbon  (CFC)  cross-sections 

-  high  altitude  weighting  functions 

-  high  spectral  resolution  (>  lcm-1)  calculations  at 

low  altitude 

-  laser/lidar  transmittances,  etc. 
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FASCODE 


PLANNED  UPGRADES: 

o  maintain  state-of-the-art  physics  and 
state-of-the-art  computational  speed 

o  new  physics  algorithms  include: 

improved  line  coupling  (C02  and  02) 
line  shape  (C02  and  H20), 

important  for  remote  sensing  applications 

o  new  radiance  algorithm  with: 

compatible  multiple  scattering, 
solar  insolation,  and  flux  capability  ; 

o  more  realistic/pratical  weighting  function  formulations 

o  AND  with  cooperation  from  DOE,  a  major  effort  to: 
implement  vectorized  coding  and  efficient/accurate 
flux  divergence  coding,  ultimately  intended  for 
parallel  processors 
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FASCODE 


USER  COMMUNITY:  There  are  over  100  new  users  of  FASCOD3P,  plus 
over  200  users  of  FASCOD2.  Users  and  uses  include: 


SSGM  (Strategic  Scene  Generation  Model)  SDIO 

TDA  (Tactical  Decision  Aids)  AF 

EOSAEL  (Army  Atmos.  Models)  Army 

SPIRITS  (Target  Surveillance  Code)  AF 

CIRRIS  data  analysis  AF 

Climate  studies  DOE 

Validation  (SHARC,  MODTRAN,  MOSART)  AF 

etc.  +200 


o  Suitable  for: 

-  all  non-solar  LOWTRAN/MODTRAN  studies 

-  recalculating  entire  systems  designs  and  algorithms 

-  simulations  relative  to  surveillance  implications 

-  implement  the  CFC  cross-sections  for  8-lOum  window  region 

-  evaluate/validate  all  systems  codes 

[any  application  that  involves  exact  or  any  broader  spectral  resolution 

under  both  equilibrium  and  non-LTE  conditions  may  be  suitable  for  FASCODE 

simulation] 
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FASCODE 


AVAILABILITY: 

FASCODE  is  available  for: 

mainframe  computers,  work  stations,  and,  potentially  PC’s. 

The  word  size  and  FORTRAN  77  coding  are  appropriately  formulated  with 
double  precision  declarations. 

FASCODE  is,  in  general,  not  being  incorporated  into  assorted  larger  codes, 
e.g.  SSGM,  EOSAEL,  SPIRITS,  TDA,  etc.,  because  of  its  required  storage  and 
running  times. 

However,  it  is  usually  maintained  "off  line"  for  validation  of 
more  pragmatic  codes. 
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M0DTRAN3,  M0DTRAN4  AND  BEYOND 
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SUMMARY 


M0DTRAIM2  /  FASC0D3P  COMMON  ELEMENTS 
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M0DTRAN2  FEATURES 
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BACKGROUNDS,  POLLUTANTS,  TARGETS,  PLUMES, 
CLIMATOLOGY,  INVERSION  ALGORITHMS,  E/0  DESIGN, 
DATA  ANALYSIS,  ETC. 


BAND  MODEL  APPROACH 


cc 

LU 

> 

o 


c/5 

LU 

I- 

< 

CC 


CD 

LU 


LU 

Q 

o 

o 

(/) 

< 


CQ  CD 


LU  ^ 
£5.  H 
CO  5 


X  Zi 
U  a. 
<  X 

LU  LU 


C/) 

LU 

a. 

< 

X 

c/5 

LU 


CD 

z 

I- 

X 


CO 

cc 

h- 

z 

o 

(J 


CD 


LL 

LU 

Q 


C/5 
LU 
H 
< 

CC 

H  CQ 

Z  Q. 
—  LU 

Z  H 
<;  c/5 

QC  LU 
l~  — J 
Q  CD 
O  Z 


CC 

LU 


H- 

15 

OQ 

CC 


CQ  h- 
5  C/5 


D 

Z 


Q 


< 

U 


LU 

> 

O  ^ 

LU  h" 
LL  C/5 

^  < 
<  CO 


LU 

CD 

< 

CC 

LU 

> 

<  f 

Z  H- 

<  CD 

Q  ^ 

<  h* 

C/5  ^ 
LU  LU 


AlIAIidUOSaV 


30 


FREQUENCY  (CM'M 


M0DTRAN2  /  FASC0D3  COMPARISON 


Ulavenumbep  C cm  - 1  > 


M0DTRAN3  UPGRADES 
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INCREASED  LAYER  RESOLUTION 


(quj)  3unss3dd 


COOLING  RATE  (K/DAY) 


SPECTRALLY  VARYING  SURFACE  ALBEDOES 
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Wavelength  (|im) 

THE  IR/EO  SYSTEMS  HANDBOOK,  Vol  1,  SOURCES  OF  RADIATION,  G.J.  ZISSIS,  Ed. 
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EXTENSION  INTO  THE  ULTRAVIOLET 
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FREQUENCY  (CM'^) 


SINGLE  SCATTER  SOLAR  RADIANCE 
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(2)  In  [4(1)  /  4(2)] 


StNGLE  SCATTER  VALIDATION 
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DISORT  MULTIPLE  SCATTERING 
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FREQUENCY  (CM'^) 


M0DTRAN3  /  FASC0D3  COMPARISON 


Ulavenumben  (cm-l) 


M0DTRAN4  AND  BEYOND 


COOLING  RATES  /  VERTICAL  FLUXES 
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-  COMBINE  COINCIDENT  &  NEARLY  COINCIDENT  LINES 

REVISE  SCATTERED  RADIANCE  CALCULATIONS 


SAMPLE  COOLING  RATE  CALCULATION 
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COOLING  RATE  (K/DAY) 


THE  ATMOSPHERE  GENERATOR  -  SAG 
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TEMPERATURE  (°K) 


SUMMARY  /  STATUS 


r 


o 

Q 

LU 

h- 

LU 

X 

H 

h- 

Q 

U 

LU 

LU 

h- 

< 

Q 

Q- 

ID 

CL 

X 

LU 

CD 

LO 

O) 

LL 

o 

CC 

LL 

CC 

(D 

O 

Z 

LU 

00 

CD 

OC 

LU 

CO 

z 

< 

LL 

Q 

LU 

G 

LU 

z 

> 

(/) 

D 

o 

CD 

oc 

z 

< 

D 

CD 

h- 

Q 

O 

z 

< 

-J 

CL 

H" 

CO 

o 

3 

_l 

LL 

o 

y 

LU 

Z 

< 

LU 

CO 

h- 

LL 

CD 

Z 

o 

o 

O 

Q 

LU  (T 

CD  < 

< 

LU 

—1 

z 

< 

LU 

<  LU 

LU 

o 

C/) 

LU 

LU  > 

CC 

Z 

UJ  CD 

z 

LU 

^  X 

z 

z 

< 

X 

1— 

< 

< 

QC 

h- 

y  cc 

CC 

CC 

H 

h- 

h- 

G 

O 

1 — 

LU 

LU 

LU 

CO  H 

?  < 

0 _ 1 

Q 

O 

G 

O 

2 

• 

• 

• 

• 

< 

< 

o 


tr 

LU 

H 

Z 

LU 

U 

< 

I- 

< 

Q 


y 

h- 

< 


-U  CN 

o  00 

CD 

O 

<  ^ 

^  ID 

CM 


O 

I- 

< 


o 


45 


QUESTIONS/PROBLEMS  WITH  MODTRAN  SHOULD 
BE  SENT  TO  GAIL  ANDERSON  OR  JIM  CHETWYND: 
FAX  (617)  377-8900 
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Hanscom  AFB,  Massachusetts 
7-8  June  1994 
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B-077-94(6-1).2 


Photon  Research  Associates,  Inc. 
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LTE/NLTE  Coupled  Code  (Pending) 


CURRENT  STATUS  OF  MOSART 
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-  Interface  with  SHARC/SAMM  Atmospheric  Generator 


ATMOSPHERIC  CHARACTERIZATION 
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Heymsfield  (Temperature  Dependent) 


GLOBAL  CLIMATOLOGY  DATA  BASES 

APRIL 


Five  Year  (1985  -  89)  Average  From  NOAA  Nimbus-7 
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BACKGROUND  REPRESENTATION 
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Galactic  and  Extra-Galactic  Radiances 


GLOBAL  TERRAIN  ALTITUDES 
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Altitudes:  -121  to  7833  Meters  at  1  Meter  Resolution 

-  Dead  Sea:  -400  m 

-  Mount  Everest:  8848  m 


SURFACE  CLASSIFICATION  MAP 
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Spatial  Resolution 


ECOSYSTEM  CLASSIFICATIONS 
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WORLD  AT  NIGHT 
(~20  MINUTE  RESOLUTION) 


PRELIMINARY  BOUNDARY  LAYER  AEROSOL  MAP 

(10  MINUTE  RESOLUTION) 
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8-077-94(6-1),  15 


OCEAN  SURFACE  TEMPERATURE 


Ref.:  S.  Levitus,  Climatological  Atlas  of  the  World  Ocean 


Photon  Research  Associates,  Inc. 
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-  User-Designed  Slit  Function 
Broad-Band  Thermai  and  Solar  Loading 
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Year-to-Year  Variations 


CODE  DOCUMENTATION 
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20%  Relative  Humidity 
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Aircraft  Reflectances  Aircraft  Reflectances 
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FUTURE  PLANS  FOR  MOSART 
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PLEXUS 


Geophysics  Directorate,  Phillips  Laboratory,  Hanscom  Air  Force  Base,  MA 

Department  of  the  Air  Force 
fax:617-377-8780,  internet  fclark@plh.af.mil 


PLEXUS 
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-  Non-interactive  hardware  independent  mode 

Proposed  upgrades 


Geophysics  Directorate  Codes 
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Dr.  Frank  O.  Clark  PL-GPOB,  HAFB 


Geophysics  Directorate  Codes 
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Geophysics  Directorate  Codes 
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Cloud  Scene  Simluation  Model 


Celestial  Background  Scene  Descriptor 
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Infrared  Astrometric  Catalog  (positions) 
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Designed  as  a  testbed 

-  determine  effective  structure  techniques 


Dr.  Frank  O.  Clark 
Geophysics  Directorate 
Phillips  Laboratory 
Hanscom  Air  Force  Base 
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weather  object  under  construction 
Unix  port  &  non-interactive  mode 
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EOSAEL  Capabilities 
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•  Support  Modules  —  PFNDAT,  AGAUS,  BITS,  CLIMAT 

•  Tactical  Decision  Aids  —  KWIK,  GRNADE,  MPLUME, 
COPTER 

•  Acoustic  Propagation  —  SCAFFIP 


TARGAC  Scenario  Specifications 


TARGAC  Selecting  Climatology  Regio 
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UVTRAIM  for  MS-WINDOWS 
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UVTRAN  for  MS-WINDOWS 


UVTRAN  for  MS-WINDOWS 
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NMMW  Gaseous  Absorption  Update 
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L/Gtailed  treatment  of  the  anisotropic,  mesospheric  Zeeman  effect  of  O2  microwave  lines 

H.  Liebe,  G.  Hufford,  and  M.  Cotton,  Propagation  modeling  of  moist  air  and  suspended 
water/ice  particles  at  frequencies  below  1000  GHz,  Proc.  AGARD  Conf.  Paper  3/1-10, 
Palma  De  Mallorca,  Spain,  May  1993. 


Verification  of  LZTRAN 
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*  H2O  and  N2  continuum  absorption  calculations. 

*  Adjustments  to  negative  extinction  coefficients. 

Sensitivity  analyses. 

*  Investigate  temp,  pressure,  water  vapor  boundary  conditions. 


Validation  of  LZTRAN 
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Examine  sensitivity  and  reasonableness  of  model  output  as  input  varies. 
Subject  matter  expert  will  review  input  and  output. 


ACOUSTICS—  SCAFFIP 
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ACOUSTICS 


H- 

i> 


+ 

iO 


II 


96 


g)soo  n  +  ( j;)3  =  s^o 


Wave  Solution 


VO 


r- 


00 


;3i 

sil 

K) 

I 


(M 

CO 


IM 


CO 


(M 


(M 

+ 


(M 

II 

ft. 

K| 

•+0 

(M 

CO 

CO 

CO 

(M 

1 

<0 

+ 

K> 

1 

CO 

CO 

a. 

(M 

c 

> 

I> 

+ 

a 

CO 

CO 

T— 1 

+ 


&H 

(M 

CO 


(M 

CO 


S7 


Complex  Ground  Impedance 
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EOSAEL  on  TECNET 
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The  Navy  Oceanic 
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THE  NAVY  OCEANIC  VERTICAL  AEROSOL  MODEL 
IN  A  FORTRAN  SUBROUTINE  FORMAT 


Stuart  G.  Gathman 


NCCOSC  RDTE  DIV  543 
53170  Woodward  Road 
San  Diego,  CA  92152 


A  FORTRAN  version  of  the  Navy  Oceanic  Vertical  Aerosol  Model,  NOVAM,  has  been 
completed  and  is  ready  to  be  interfaced  with  other  types  of  transmission  models.  The  model  is 
available  in  the  form  of  a  NCCOSC  RDT&E  DIV  technical  report  #1634  which  contains:  (a)  a 
written  description  of  the  model,  (2)  the  printed  code  of  the  subroutine,  (3)  the  printed  code  of 
the  associated  subroutines  and  functions,  and  (4)  a  sample  driver  which  can  be  used  for  batch 
processing  of  data.  The  model  requires  meteorological  data  in  the  form  of  both  surface  and 
atmospheric  sounding  data  files.  NOVAM  predicts  the  venical  distribution  of  aerosol  in  the  first 
6  km  above  the  ocean.  Outputs  of  the  model  include  optical  properties  in  the  wavelength  band 
from  0.2  to  40  micrometers  at  any  altitude.  These  properties  include  the  volume  extinction  and 
absorption  coefficient,  the  relative  humidity,  as  well  as  parameters  describing  the  aerosol  size 
distribution  in  terms  of  dN/dr  at  that  altitude.  This  paper  will  describe  the  model  and  some 
sample  applications. 
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Only  available  meteorology  data  to  be  used. 
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earlier  methods  not  really  satisfactory  for 
LOWTRAN/MODTRAN 

use  of  visibility  to  determine  a.m.p. 
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Surface  Observation  file 
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Temperature 


Mixing  Ratio  profiles  recognized  by  NOVAM 
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EXECUTIVE  SUMMARY 


OBJECTIVE 

The  object  was  the  development  of  the  FORTRAN  version  of  the  Navy  oceanic  vertical 
model  (NOVAM).  The  model  predicts  the  vertical  distribution  of  aerosol  in  the  first  6000  meters 

above  the  ocean. 

RESULTS 

The  NOVAM  was  developed  from  extensive  marine  aerosol  studies  from  different 
laboratories.  The  climax  was  a  multimenu-driven  interactive  program  that  allows  mouse 
selection  of  menu  items  needed  for  the  calculation.  It  includes  graphics  and  editing  capabilities 
useful  to  the  researcher.  When  used  with  an  appropriate  method  to  determine  the  profile 
parameters,  it  could  be  used  in  a  fully  automatic  mode  in  which  the  program  could  access  sets  of 
data  and  produce  an  analysis  of  the  data  sets  in  an  unattended  background  mode.  The  final 
NOVAM  code  is  intended,  however,  to  be  used  in  conjunction  with  a  LOWTRAN/MODTRAN 
program  and  to  supply  the  electro-optical  (EO)  propagation  characteristics  to  the  calling 
program  that  are  produced  by  the  unique  aerosol  found  in  the  marine  atmosphere.  It  is  written  in 
FORTRAN  so  it  can  be  integrated  into  LOWTRAN/MODTRAN  codes  to  improve  model 
performance  in  marine  environment. 

RECOMMENDATIONS 

The  model  has  several  shortcomings  that  will  be  addressed  in  future  modifications.  The 
region  of  applicability  leaves  two  areas  not  covered  well  by  the  model.  First,  higher  altitudes, 
various  models  developed  by  the  U.  S.  Air  Force  and  included  in  the  LOWTRAN/MODTR.AN 
code  will  be  more  accurate.  Second,  propagation  paths  that  graze  the  sea  surface  or  pass  through 
the  region  within  a  meter  or  so  of  the  sea  surface  are  not  adequately  covered  by  NOVAM.  This 
problem  is  being  remedied  by  a  large-scale  experiment  off  the  Dutch  coast  sponsored  by  NATO. 
Another  shortcoming  is  the  limited  types  of  weather  situations  in  which  it  is  applicable. 
Advances  in  these  areas  await  development  of  models  from  the  basic  research  community  in  t  e 
future  Another  area  of  concern  is  the  use  of  the  model  in  close-in  coastal  areas.  Compensation 
was  introduced,  but  experience  has  shown  this  is  one  of  the  weakest  parts  of  the  model.  It  is  the 
author’s  opinion  that  a  special  coastal  aerosol  model  needs  to  be  developed  that  will  adequate  y 
take  into  account  local  sources  of  aerosol. 


CONCLUSIONS 


We  have  presented  a  model  for  describing  the  EO  propenies  of  the  unique  marine  aerosol 
found  in  the  regions  from  shipboard  height  to  above  several  kilometers  in  altitude.  The  model 
has  been  written  as  a  self-contained  FORTRAN  subroutine  so  it  could  be  incorporated  into  larger 
scale  models  such  as  the  LOWTRAN  and  MODTRAN  codes.  The  model  needs  information  on 
the  meteorological  sounding  at  the  site  where  the  calculation  is  made  as  well  as  information  on 
certain  meteorological  parameters  near  the  surface  of  the  sea.  The  model  has  certain  shortcom¬ 
ings  that  need  to  be  addressed  in  future  modifications.  First  of  all.  the  region  of  the  applicability 
is  from  shipboard  level  (about  5  to  10  meters)  to  regions  abo\’c  the  lower  troposphere  where  oth¬ 
er  aerosol  models  will  be  more  appropriate.  This  leaves  two  areas  that  are  not  covered  well  by 
the  model.  At  the  higher  altitudes,  various  models  developed  by  the  U.  S.  Air  Force  and  included 
in  the  LOWTRAN/MODTRAN  codes  will  be  more  accurate.  On  the  other  extreme,  an  important 
propagation  path  that  grazes  the  sea  surface  or  passes  through  the  region  within  a  meter  or  so  of 
the  sea  surface  is  not  adequately  covered  by  NOVAM.  This  is  because  NOVAM  is  in  part  an 
empirical  model  and  based  on  measurements  in  the  real  world.  Current  interest  from  shipboard 
level  on  down  currently  lacks  observation  data  because  of  the  difficulty  in  obtaining  them.  This 
region  will  be  especially  important  to  IR  propagation  during  rough  weather  and  high  seas  where 
many  marine-generated  drops  and  droplets  are  suspended  in  these  lower  levels  of  the 
atmosphere.  This  problem  is  currently  being  remedied  by  a  large-scale  experiment  called  the 
Marine  Aerosol  Properties  and  Imager  Performance  (MAPTIP)  trial  off  the  Dutch  coast 
sponsored  by  NATO.  The  results  of  this  experiment  will  contribute  to  the  development  of  an 
advanced  Navy  aerosol  model  (ANAM)  currently  under  development.  These  results  can  be 
added  into  the  modular  format  of  NOVAM  to  increase  its  regions  of  application. 

Another  shortcoming  in  NOVAM  is  its  somewhat  limited  types  of  weather  situations  in 
which  it  is  applicable.  An  earlier  version  of  NOVAM  included  the  region  just  below  stratus 
clouds,  but  because  this  model  had  a  limited  band  of  wavelength  validity  and  required  inputs  that 
are  really  incompatible  with  a  self-contained  model,  this  submodel  was  dropped  from  the  current 
model.  Advances  in  these  areas  await  the  development  of  models  from  the  basic  research  com¬ 
munity  sometime  in  the  future. 

An  area  of  concern  in  the  application  of  models  such  as  NAM  and  NOVAM  is  the  use  of  the 
models  in  the  close-in  coastal  areas.  As  these  models  were  developed  for  the  open  ocean  region 
far  away  from  the  land  influences,  error  would  be  expected  when  unusual  sources  of  aerosol  are 
sent  into  the  atmosphere  by  man-made  sources.  The  amp  concept  was  introduced  into  NAM  and 
NOVAM  to  compensate  for  these  problems,  but  experience  has  shown  that  this  is  one  of  the 
weakest  parts  of  the  models.  It  is  the  author’s  opinion  that  a  special  coastal  aerosol  model  needs 
to  be  developed  that  will  adequately  take  into  account  local  sources  of  aerosol. 


THE  NATURE  OF  DESERT  AEROSOLS:  SUMMARY  OF 
CHINA  LAKE  STUDIES  TO  DATE 


P.L.  Walker 

Physics  Department,  Code  PH 
Naval  Postgraduate  School 
Monterey,  CA  93943-5117 
walker(2)physics.nps.navy.mil 

L.A.  Mathews,  ret. 

Research  Department 

Naval  Air  Warfare  Center  -  Weapons  Division 
China  Lake,  CA  93555 

ABSTRACT 

The  Naval  Air  Warfare  Center  has  been  conducting  desert  aerosol  characterization  measurements  at 
China  Lake  at  ground  level  in  1987;  at  15,000  feet  in  connection  with  Long  Jump  in  1988;  at  ground 
level  in  the  western  Mojave  including  Edwards  AFB  and  China  Lake  in  1990;  and  will  be  conducting 
further  measurements  in  the  western  Mojave  in  the  summer  of  1994.  This  presentation  will  present 
an  over  view  of  the  analysis  of  data  to  date.  The  properties  of  desert  aerosol  that  are  appearing  are 
some  what  different  than  those  proposed  by  Longtin  and  Shettle  in  that  the  accumulation  mode  is 
dominated  by  the  presence  of  large  numbers  of  organic  carbon  particles  while  dust  is  characterized 
by  the  presence  of  clay  particles  in  the  dust  size  distribution  mode  with  few  quartz  particles  being 
present  and  no  quartz  dominated  blowing  sand  mode. .  We  will  conpare  particle  size  and  composition 
at  15,000  feet  and  ground  level,  further  micrographic  studies  of  the  composition  of  high  altitude 
aerosols  and  the  effect  of  wind  speed  on  aerosol  characteristics  at  ground  level. 
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MEASUREMENTS 
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•  Stage  9:  0.1  to  0.4  jim 


A  Comparison  of  Ground  Aerosol  Sizes  and  Concenlralions  with  Those  Near  15,500  ft. 

Dust  Mode 


@15,500  ft. 


@  Ground 


@15,500  ft. 


@Ground 


No  =  0.6  particles/cc  number  density 

Dn  =  0.285  microns  mode  diameter 

Log(o)  =  0.7  mode  width 

No  =  5  particles/cc 
Dn  =  0.35  microns 
Log(a)  =  0.9 


Accumulation  Mode 

No  ~  5  X  lOVcc 
Dn  =  .07  p 
Log(o)  -  0.28 

No  ~  2  X  10’ /cc 
Dn  =  s  O.lp 
Log(o)  ^  0.47 
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Accumulation  Mode  Particle  Concentration  vs  Altitude  for  Flight  15. 


Particles/cc 


Particle  Composition  at  15,000  ft 

•  Calcium  and  phosphorus  particles  having  featureless  surface 
characteristics  were  found  in  all  size  ranges.  One  particle 
was  associated  with  chlorine,  sulfur  and  possibly  clay  in  the 
0.1  to  0.4p  size  range  during  stagnant  atmospheric 
conditions. 

•  Clay  was  found  in  all  size  ranges  and  always  in  association 
with  sulfur  for  stagnant  conditions,  but  only  for  the  0.4  to 
l.Op  range  during  storm  conditions. 

•  Sulfur  is  not  restricted  to  the  accumulation  mode.  Sulfur 
composes  an  increasing  proportion  of  the  alumino¬ 
silicate/sulfur  mass  for  decreasing  particle  size.  This 
phenomena  is  consistent  with  the  concept  of  sulfuric  acid 
absorption  into  alumino-silicate  surface  layers. 

•  Chlorine  was  also  associated  with  alumino-silicates  in  all  size 
ranges  for  stagnant  conditions  but  was  not  sought  for  storm 
conditions.  This  association  is  consistent  with  accumulation 
of  chlorine  from  sea  salt  under  stagnant  atmospheric 
conditions. 

•  Chlorine  and  sulfur  occurred  in  a  fixed  ratio  except  for  one 
instance. 

•  Gypsum  was  found  in  size  ranges  greater  than  one  micron  for 
stagnant  conditions. 

•  Sodium  sulfate  was  found  in  the  0.4  to  l.Op  range  during 
storm  conditions. 
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Smog  and  Dust  in  the  Mojave  Desert 


Ground  level  atmospheric  extinction  from  0.5  to  12  microns  was  determined 
as  a  function  of  date  and  time  of  day  by  a  combination  of  direct  visibihty 
measurements  and  Mie  calculations. 

Simultaneous  Measurements  were  made  using 

Nephelometer 

Telephotometer 

Aerosol 

instruments  in  the  Indian  Wells  Valley  from  April  through  July  1987. 

Spectral  Extinction  was  calculated  using  measured  aerosol  size  distributions 
and  composition  estimated  from  the  RESOLVE  Report. 

Results 

Unexpectedly  high  extinction  in  the  8  -  12  micron  range  agrees  with 
anecdotal  information:  In  side  by  side  operation  at  the  Naval  Weapons 
Center  detectors  operating  in  the  8  -  12  do  not  perform  as  well  as  those 
operating  in  the  3  -  5  micron  range. 

The  situation  may  be  worse  than  what  we  have  calculated: 

Mid-day  (when  most  of  our  data  were  taken)  accumulation  mode  aerosol 
concentrations  were  only  60%  of  the  twenty-four  hour  average  and 
sometimes  only  a  quarter  of  the  midnight  concentrations. 

Measured  organic  aerosol  concentration  may  be  too  high  by  a  factor  of  two 
due  to  2X4  filter  contamination.  That  would  imply  that  (NH4)2S04  induced 
extinction  could  be  twice  what  we  have  calculated. 


Time  of  Day  (PST) 


FIGURE  21.  Daytime  Dust  Mode  Concentrations.  The  middle  points  were 
averaged  over  three  months  starting  in  May  1987. 
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Table  8  .  Monthly  Average  Log-Normal  Fit  Parameters. 


Organic  Fine  Mass 

Month  Soot  Sulfate  Nitrate  Carbon  Dust  (mcrg/m''3) 


REPRODUCED  AT  GOVERNMENT  EXPENSE 
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Material 


Montmorillonite  Clay 
Kaolinite  Clay 
lllite  Clay 
Quartz 

Sodium  Nitrate 
Calcite 

Potassium  Nitrate 
Dolomite 
Hematite 
Halite 


%  of  Dust 

Density 

3.57 

2.5 

5.32 

2.64 

25.3 

2.75 

22.38 

2.65 

23.45 

2.26 

2.41 

2.71 

2.06 

2.11 

5.79 

2.86 

5.71 

5.24 

4.01 

2.17 

Average  dust  density  =  2.61  gm/cc 
Table  17.  China  Lake  Dust  Using  Pye  Model. 


Material  %  of  Dust 


Montmorillonite  Clay 
Kaolinite  Clay 
lllite  Clay 
Quartz 

Sodium  Nitrate 
Calcite 

Potassium  Nitrate 
Dolomite 
Hematite 
Halite 

Average  dust  density  =  2.55  gm/cc 

Table  18.  China  Lake  Dust  Using  Hoidale  Model. 


20.3 

10.1 

10.6 

15.5 

21.1 

3.6 

5.6 
4.2 
5.4 
3.8 
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I>ust  Mass  (mcgm/te**3) 


Figure  ,  Dust  Mass  versus  Windspeed. 
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Dust  Mode  Diameter,  microns 


Figure  .  Dust  Mode  Diameter  versus  Windspeed. 
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Log-Width  of  the  dust  Mode 
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Mass  Distribution  from  beginning  to  end  of  Storm 
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Windspeed  and  Dust  Loading  at  Edward 


ve  Hour  Averaged  Wind  Speed 


Silicon  to  Aluminum  Mass  Ratio  vs  Wind  Speed.  The  authors  do  not 
know  why  the  silicon  to  aluminum  content  decreases  with  wind  speed. 


Future  Work 


Investigate  organic  aerosol  problem  by  comparing  mass  captured  on 
filters  to  mass  calculated  from  size  distribution.  Use  a  denuder 
(pre-filter)  to  remove  organic  vapors. 


Use  LIDAR  to  get  vertical  aerosol  profiles  to  15,000  to  20,000  ft. 


Make  greater  effort  to  get  all  the  instruments  working  all  the  time  in 
order  to  obtain  comparison  and  rare  event  data. 


Use  impactor  and  x-ray  elemental  mapping  more  extensively  than 
here-to-fore  in  order  to  deduce  particle  composition  statistics. 


Relate  aerosol  and  meteorological  data  to  meteorological  codes  that 
will  run  on  a  work  station.  Codes  that  do  not  appear  to  be  suitable 
are  NOGAPS,  CO  AMPS  and  NORAPS  since  they  take  hours  to  run 
on  a  Cray  C-90. 
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Application  of  the  MODTRAN2  code  to  the  modeling  of 
silicate  dust  clouds. 

Abstract 

The  MODTRAN2  code  has  been  used  to  model  the  thermal 
infrared  radiance  from  a  dust  cloud  formed  during  an  explosives  test. 
High  signal-to-noise  ratio  infrared  spectra  (3  -  14  microns,  XI AX  =  0.1 
microns)  of  the  dust  cloud  were  measured  and  they  showed  the 
presence  of  silicate  particles.  In  an  attempt  to  model  the  thermal 
emission  of  the  dust  cloud  and  reproduce  the  silicate  feature,  a 
silicate  dust  cloud  was  generated,  and  MODTRAN2  was  used  to 
calculate  radiance  spectra  between  6  and  16  microns.  A  description 
of  the  analysis,  comparison  with  experimental  data,  and  the 
problems  encountered  will  be  presented. 


1.  Observations 

An  explosives  test  in  U.S.  southwest  desert  lofted  a  large 
amount  of  dust  into  the  atmosphere  to  altitudes  of  several  thousand 
meters  above  ground  level  (AGL).  High  signal  to  noise  spectra  of  the 
resulting  dust  cloud  were  measured  in  the  3-14  micron  range  using  a 
nonscanning  array  spectrograph  (Hackwell  et  al.  1990).  Data  were 
taken  during  the  initiation,  growth  and  dispersion  of  the  dust  cloud 
(Lynch  et  al.  1994).  In  addition,  spectra  of  the  clear  sky  and  that  of 
a  cumuloform  water  cloud  were  also  collected.  The  spectra  of  the 
dust  cloud  showed  clear  evidence  of  the  silicate  SixOy  vibrational 
emission  band  near  10  microns  (Figures  1,  2),  Laboratory  analysis  of 
soil  samples  collected  from  ground  zero  before  and  after  the  test 
confirmed  the  identification. 

Owing  to  the  limited  slant  angles,  column  densities,  composition 
and  other  atmospheric  variables,  we  decided  to  model  the  dust  cloud 
radiance  using  a  general  purpose  atmospheric  radiance  code.  The 
initial  objective  was  to  calculate  the  infrared  spectra  expected  from  a 
model  silicate  cloud  and  compare  this  to  the  observed  data.  The 
MODTRAN2  model  was  chosen  based  on  available  expertise  in 
running  the  code,  adequate  spectral  resolution  to  match  the 
instrument  used  for  the  observations,  and  an  existing  dust  model 
built  into  the  code. 
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The  very  existence  of  the  silicate  emission  feature  gives  some 
indication  of  the  particle  size.  Optically  thick  particles  would  not 
show  a  silicate  feature.  For  typical  silicates  for  which  the  imaginary 
parts  of  the  indices  of  refraction  around  10  microns  are  near  unity 
(highly  absorbing),  the  maximum  particle  size  that  would  show 
silicate  emission  is  of  order  10  microns  in  radius.  With  this  in  mind, 
we  used  the  indices  of  refraction  for  various  silicate  minerals 
(including  MODTRAN2’s  volcanic  aerosols),  computed  the  absorption 
and  extinction  efficiencies  Qabs  and  Qext  using  Mie  theory,  then 
defined  a  dust  cloud  at  a  height  and  size  consistent  with  the  cloud 
produced  in  the  test.  This  cloud  was  approximately  500  meters  thick, 
approximately  3  km  AGL  and  observed  from  the  ground  with  a 
zenith  angle  of  about  60  degrees.  The  extinction  coefficients  were 
computed  and  used  as  inputs  to  MODTRAN2. 


3.  Aerosols  and  MODTRAN2 

In  operation,  the  MODTRAN2  model  will  define  a  model 
atmosphere  layering  based  on  predefined  altitudes  for  four  different 
aerosol  regions.  A  first  attempt  at  defining  a  model  atmosphere  was 
made  by  allowing  MODTRAN2  to  generate  its  own  layering.  For  this 
observing  geometry,  and  given  our  plans  to  change  the  defined 
aerosol  layer  contents,  the  default  layering  was  inadequate.  While 
MODTRAN2  performs  some  scaling  of  the  aerosol  profiles  based  on 
the  actual  starting  altitude  for  the  observation  (through  the  GNDALT 
input),  the  boundaries  of  the  four  built-in  aerosol  layers  are  fixed  at 
0  to  2  km,  2  to  10  km,  10  to  30  km,  and  30  to  100  km.  For  this 
viewing  site,  the  altitude  for  our  test  was  close  to  the  first  aerosol 
layer  boundary  at  two  kilometers.  Because  of  this,  the  MODTRAN2 
atmosphere  allowed  only  three  layers  for  the  boundary  layer  aerosol 
before  the  second  predefined  aerosol  layer  boundary  was  reached. 

Using  the  MODTRAN2  generated  layering  as  a  guide,  a  model 
atmosphere  was  created  using  100  meter  laying  through  4 
kilometers,  with  the  remaining  layers  spread  upwards  to  100 
kilometers.  This  was  the  only  way  to  redefine  the  aerosol  layer 
boundaries  in  the  MODTRAN2  inputs.  One  of  the  difficulties  in  this 
process  occuixed  because  MODTRAN2  limits  the  model  atmosphere  to 
only  34  layers.  In  this  case  the  layer  limitations  required  that  the 
stratospheric  ozone  layering  be  made  coarser  to  allow  for  better 
modelling  of  the  boundry  layer  aerosols. 
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The  extinction  and  scattering  coefficients  of  the  dust  were 
calculated  using  the  Mie  theory  code  developed  by  Bohren  and 
Huffnnan  (1983).  This  calculation  assumes  that  the  particles  are 
spherical  in  shape.  While  the  assumption  of  sphericity  is  valid  for 
liquid  aerosol  particles,  the  variety  of  shapes  in  dust  particles  can 
give  rise  to  errors.  Although  this  assumption  may  yield  answers  that 
are  not  entirely  correct,  the  uncertainties  in  the  optical  constants  and 
size  distributions  are  certainly  more  significant.  The  values  of  Qext 
and  Qscat  were  calculated  for  each  of  the  47  wavelength  values 
required  by  the  MODTRAN2  code.  The  optical  constants  for  the 
silicate  material  were  obtained  from  Palik  (1982),  which  compiles 
the  data  from  numerous  sources. 

In  addition  to  using  these  calculated  values,  the  extinction  and 
absorption  coefficients  for  the  volcanic  aerosols  available  in 
MODTRAN2  were  used  to  define  the  dust  cloud.  The  structure  in  the 
’aged  volcanic'  absorption  suggests  that  this  aerosol  contains  some 
silicates  (Volz  1973,  Shettle  and  Volz,  1976). 

The  scattering  coefficients  were  then  included  into  the  second 
of  the  four  aerosol  regions  in  the  model  atmosphere  described  above. 
Only  the  aerosol  species  of  interest  was  changed  when  the  dust  cloud 
was  loaded  into  the  model  atmosphere;  the  other  aerosol  species  and 
the  atmospheric  layering  remained  unchanged. 


4.  Comparison  with  Observations. 

Figure  3  shows  a  comparison  between  the  observed  spectrum 
and  that  calculated  by  the  default  layering  generated  by  MODTRAN. 
The  tropical  atmosphere  profile  was  required  to  match  the  observed 
radiance,  indicating  a  large  amount  of  water  vapor  was  present  in 
the  atmosphere  even  though  the  observing  site  was  in  a  desert.  The 
deficit  between  the  calculated  and  observed  radiances  is  probably 
due  to  a  hygroscopic  aerosol  haze  in  the  atmosphere,  and  could  be 
modelled  by  using  a  rural  aerosol  model.  Further  investigation  of 
this  will  be  done  using  the  sonde  data  in  a  follow-on  study. 

Figure  4  shows  a  comparison  between  the  observed  radiance 
and  a  MODTRAN2  calculated  spectra  using  the  'aged  volcanic’  aerosol 
absorption  and  extinction  coefficients  that  are  built  into  the 
MODTRAN2  code.  The  observed  spectrum  used  for  this  comparison 
was  taken  thirty  minutes  after  the  initial  cloud  was  formed,  so  the 


146 


dust  particles  should  be  in  thermal  equilibrium  with  the  surrounding 
atmosphere.  The  agreement  between  the  two  is  reasonable,  showing 
the  expected  silicate  emission  in  the  8.5  to  9.5  micron  band. 
Unfortunately,  since  the  actual  composition  and  particle  sizes  of  the 
'aged  volcanic'  material  were  not  known  to  us  at  the  time  of  this 
study,  a  further  comparison  could  not  be  made. 

In  Figure  5  the  radiance  from  a  cloud  composed  only  of 
amorphous  Si02  particles  of  5  micron  radius  is  compared  with  the 
observed  cloud  data.  The  differences  in  these  spectra  between  8.5 
and  9.5  microns  clearly  indicate  that  the  structure  of  the  amorphous 
silicate  emission  is  not  adequate  to  model  this  cloud.  This 
comparison,  combined  with  laboratory  analysis  of  dust  from  the  site, 
suggested  the  use  of  a  crystalline  silicate  material  to  better  model 
the  spectral  structure. 


5.  The  47  wavelength  problem. 

When  we  began  to  include  a  crystalline  silicate  material, 
quartz,  into  a  dust  cloud  model,  a  significant  problem  was 
encountered.  Figures  6  and  7  show  the  scattering  coefficients 
calculated  for  input  into  MODTRAN2  for  a  5  micron  radius  cloud 
composed  of  quartz,  using  the  optical  constants  from  Palik  (1985). 

The  full  resolution  values  are  drawn  in  the  solid  lines,  while  the 
values  which  MODTRAN2  uses  are  shown  as  the  dashed  line.  Note 
that  significant  spectral  structure  is  lost,  particularly  near  the  9.5 
micron  absorption,  as  well  as  near  the  12.5  micron  absorption  bands. 

Apart  from  selecting  which  optical  constants  to  use  for  the  dust 
cloud  model,  the  fixing  of  the  wavelengths  at  which  the  optical 
constants  can  be  input  to  the  MODTRAN2  code  poses  a  serious 
limitation  to  accurately  modelling  the  emission  from  the  dust  cloud. 
While  the  model  does  allow  for  an  input  wavelength  from  card  2D2, 
these  values  are  discarded  since  they  must  conform  to  the  required 
47  input  wavelengths.  While  the  choice  of  these  47  wavelengths 
does  perhaps  make  some  sense  in  the  context  of  an  overall 
atmosphere  (e.g.  interference  from  the  ozone  band  near  9.5  microns), 
this  imposes  a  limitation  on  the  types  of  materials  that  can  be 
modelled.  This  limitation  does  not  render  the  model  useless,  since 
these  wavelengths  do  work  well  for  some  materials.  For  example 
figures  8  and  9  show  a  comparison  of  the  amorphous  Si02  optical 
constants  that  can  be  input  to  MODTRAN2  with  the  actual  values.  For 
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this  material,  the  fixed  wavelengths  are  a  reasonable  approximation 
to  the  real  values. 


6.  Discussion 

In  attempting  to  define  our  own  aerosol  cloud,  a  number  of 
difficulties  were  encountered  in  using  MODTRAN2  code.  While  the 
LOWTRAN  manuals  provided  some  guidance  for  constructing  a  model 
atmosphere,  the  process  of  introducing  a  user  defined  aerosol 
brought  out  several  questions.  Among  these,  there  is  a  question  as  to 
why  there  are  four  regions  allowed  for  the  user  defined  aerosols,  yet 
there  are  five  total  aerosol  regions  which  are  defined  in  the  code. 

The  use  of  this  fifth  aerosol  region,  and  therefore  its  effect  on  the 
data,  is  still  unclear. 

A  second  issue  is  the  composition  and  particle  size  distributions 
of  the  built-in  dust  aerosols.  Volz  (1972)  and  Shettle  and  Volz 
(1976)  discuss  the  absorption  and  transmission  of  field  gathered  dust 
particles  without  a  discussion  of  the  composition  of  dust  itself. 

As  discussed  above,  the  restriction  of  wavelengths  at  which  the 
scattering  parameters  can  be  input  limits  the  aerosol  materials  that 
can  be  modelled.  To  explore  this  further,  we  have  modified  the 
MODTRAN2  code  to  remove  this  wavelength  restriction.  However  the 
restriction  to  47  wavelength  values  has  not  yet  been  removed,  and 
makes  the  selection  of  the  appropriate  wavelengths  for  each  aerosol 
material  very  important.  The  results  of  this  investigation  will  be 
discussed  in  a  future  study.  It  is  interesting  to  note  that  other 
investigators  have  also  encountered  this  difficulty.  Hammer  has 
made  similar  modifications  to  the  LOWTRAN  code  in  connection  with 
his  work  on  cirrus  clouds  (Hammer  et  al.  1994). 


7.  Conclusions  and  directions  for  further  work 

The  MODTRAN2  model  works  reasonably  well  for  modelling 
the  infrared  spectra  from  a  silicate  cloud  in  a  model  atmosphere. 


The  user  defined  aerosol  layer  capability  of  the  MODTRAN2 
model  works  well  for  materials  whose  spectral  structure  is  smooth  so 


that  the  scattering  parameters  are  well  sampled  by  the  available 
input  values. 

We  recommend  the  following  modifications  of  MODTRAN2  be 
made.  The  existing  restriction  of  47  fixed  wavelengths  at  which  the 
scattering  parameters  can  be  input  should  be  removed.  This 
restriction  limits  the  aerosol  materials  that  can  be  modelled  using  the 
MODTRAN2  code. 
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THE  AEROSPACE 
CORPORATION 


_ - — - — - 1 

QNERA 

Planche;  1 

ONERA 

AirbornG  iriGasurGmGnts  of  cloud  radiation 

and 

comparison  with  theory 

C.Malherbe,  P.Simoneau,  P.Michon 
A.Boischot,  G.Durand,  J.Deschamp,  G.Gregoire 


1/  Introduction 

2/  Airborne  measurements 

3/  Comparisons  with  LOWTRAN7 

4/  Improvement  of  the  model  ! 
5/  Results 

6/  Conclusion 
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Planche: 


INTRODUCTION 


Requirements  at  ON  ERA 


Modelisation  of  the  infrared  atmospheric  and  terrestrial  backgrounds 
in  presence  of  clouds: 

-  fluctuations  of  the  atmospheric  background 

-  fluctuations  of  satellites  terrestrial  pictures 


Tools 


-  airborne  measurements 

-  computer  models 


164 


CLOUD  CAMPAIGN  ONERA 


ONERA 


Planche:  3 


ONERA 


05  cC  cC 

N  N  N 
C/)  CO  03 

000 


SICAP  Spectra  (Spectrometre  Infrarouge  Cryogenique  AeroPorte) 


Planche:  4 

Airborne  measurements 


Spectrometre  Infrarouge  Cryogenique  AeroPorte  (SICAP) 
(airborne  cryogenic  IR  spectrometer) 


IR  Spectrometer:  -  circular  variable  filter 

-  band:  1 .5  -  5.5  jxm 

-  resolution:  AX/X  =  2% 


>  Cryogenic 


low  thermal  noise 


Airborne 


C.E.V.'s  caravelle  116 


•3.66 


Planche:  SICAP 


Spectrometre  Infrarouge 
Cryogenique  AeroPorte 

(SICAP) 


088 
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Longitude  ~  deg 


Cloud  campaign:  2  july  1992  /  SZA  =  4B  deg. 


13 


iOSZA  =  107  (kg.  /  Rel.  Mm.  =  132  (kg.  /  Scat,  ingle  =  132  (kg. 


'4 


Planche:  7 


Discussion 


Comparisons  measurements  with  LOWTRAN7 

— ►  approximative  method  for  the  computation 

of  the  multiple  scattering  ( two  stream  method  ) 


Improvements 

— ►  more  accurate  method  (  Discrete  Ordinate  Method  ) 
— ►  microphysical  parameters 


NUALUM 
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OMERA 


Planche:  8 


Cloud  radiation  pattern 


n*F(i 


Ltm(A,,Qcp) 

T 


Iot(X,,Q,(pb) 


V 


lcloud(A,,Q(p) 


ltot(A,,Qi,cp))  =  jLlo  n 
+  Tco(A) 


Fo  R(A,,6b,Cp),€b,(pb)  Xco(A)+  lcloud(A/,Q,C[)))  Tco(A,)+  latm(A<,6b,<p)) 

I  R(A,eb,Cp),6,C^  latm(?l,6,cjW0Clcp 
•^e,cp 


(io  =COS0b 

lcioud(A,,0,cpD)  =  S{X,03)  *  B(A,T) 


R(^,e,cp,Q(p'),  eix,ei,  itfo,  b(a,t),  iatm(>.,e,() ,  Tco(?i) 


7 
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QNERA 

iC 

Planche:  9 

QNERA 

Organigram  (i) 


1/  Mie  scattering  code  (J.V.Dave) 

kext(i,?i)  =J*dn(i,r)/dr  C5^xt(A)  dr 

kabs(i,^)  =J*dn(i,r)/dr  cfebs(A)  dr 

g(i,5i)=1/2  I  P(i,cos0)cos0dcos9 
•'-1 

=  ke)(t(i,A)  Az 


for  each  layer "  i "  in  the  cloud  ( inhomogeneous  cloud) 
and  n(r)  given  by  Diem's  distribution  with:  -  r  =  5  |um 

-  eqlwc  =  0.6  g/cm^ 


2/  DISORT  (  K.Stamnes  and  Collaborators) 

(  discrete  ordinate  radiative  transfert  ) 

kext(i,A-) 

kabs(i,A,) 


DISORT 


32  stream 


R(X,0,cp,$(p’ 

e(X,e) 
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ONERA. 

17 

Planche;  1 0 

ONERA 

Organigram  (2) 


3/  LOWTRAN7  (PLYGL) 


Atmospheric  data 
( radiosoundings) 

P(z),T(z),  [H2O] 
Ztop,  Zbottom 


LOWTRAN7 


4/  NUALUM  (ONERA) 

{  nuage  luminance  ) 


T„(X) 

nFo 

I  atm 

R(X.,0,(p,Qcp’ 

e(X,e) 


NUALUM 


nFo 
►  T»(X) 

I  atm  (X,,e,(p)) 


ltot(X/,Q,(p)) 


with  J*  R(?i,0,cp,Qcp’)  latm(}i,0,cpl)  d0dcp'  « 1  lFo*R(X.,0,cp,Qcp') 
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Z.sr.microns 


microns 


LOSU  -  I  OS  de§.  /  M.  him.  =  42  deg.  /  5co(,  Angle  =  70  deg, 


fr/cm>2.  snmicrons 


Planche:  1 1 


CONCLUSION 

Fair  agreement  between  measurements  and  models  (  LOWTRAN7, 
and  NUALUM) 

Improvement  of  the  multiple  scattering  model  by  DOM 

PERSPECTIVES 

Constitution  of  a  data  bank  from  DISORT  for  different  types  of  clouds 
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UV  PROPAGATION  PHENOMENOLOGY 
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w  Ozone  Low  Ozone  High  Ozone  High  Ozone 

w  Scatter  High  Scatter  Low  Scatter  High  Scatter 


MWR  MODEL 
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TPART  MONTE  CARLO  (RADIATION  RESEARCH  ASSOCIATES) 
ZACHOR  RECURSION  FORMULA 
RIEWE  &  GREEN  MONTE  CARLO 
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SCATTERING  ANGLE 


Inclusion  of  Hccurate  Multiple  Scattering  in 

MODTRRN 
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calculations  of  atmospheric  transmission  and  prediction 

accuracy  of  retrieuals  of  remoteiy-sensed  atmospheric  properties  that  rely 
on  the  interpretation  of  measured  radiances  (ground-based  or  from  space) 
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*  extensive  testing  is  required  and  must  be  carefully  executed.  This  is  a  time- 
consuming  undertaking. 
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uery  efficient.  Therefore  it  luould  be  of  great  interest  to 
find  out  under  uihat  conditions  it  is  ualid. 


B.  Multiple  Scattering  in  Spherical  Geometry 
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B.  For  isotropic  scattering  there  is  nothing  'driuing'  azimuth-dependence.  So  the  solution  iiiill  be 
azimuth-independent  in  this  case. 


.  Multiple  Uieujing  Directions 
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radiance  computation. 
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path  consists  of  sections  with  uniform  properties  the  integral 
may  be  expressed  as  the  exponential  of  a  double  sum. 
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Eye-Safe  Laser  Ranger  Options 
Available  For  Airborne  Sensor  Program 
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Laser  Atmospheric  Extinction  vs  Altitude 

FASCOD3P/HITRAN92:  MLS,  Bkg  Strat/Mod  Vole,  Rural-23km 
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One-Way  Laser  Atmospheric  Transmission 

500  km  range  to  target  from  35  kft  Aircraft 
MLS,  Bkg  Strat/Mod  Vole,  Rural-23km  Atmospheric  Model 


0>eOI^<OU)^CO<Mir-0 

■  ■•■■■••a 

ooooooooo 


aoueuiuisueji  Ae/w-^uo 


284 


Target  Altitude  (km) 


Transmittance  Using  Standard 
LOWTRAN  Stratospheric  Aerosol  Models 
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Why  SAGE  II 

(Stratospheric  Aerosol  and  Gas  Experiment) 
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Hampton  (LaRC)  Lidar  Parameters 
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Mount  Pinatubo  6/12/91 
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Middle  East  SAGE  profiles 
1 984  to  1 992 
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ddle  East 


Aerosol  Extinction  Coefficients  (km"’ 


SAGE  1985  Middle  East 
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LaRC  Lidar  Scattering  Ratios 
8/3/91  to  2/19/93 
Pinatubo  Erupted  in  June  91 
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ref.  M.  Osborn,  private  communications  1993 


SAGE  Extinction  Profiies  over  Hampton 
from  5/10/91  to  12/31/92 
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The  number  above  each  figure  gives  the  date.  The  first  two  numbers  specify  the 
year,  the  second  two  the  month,  and  the  iast  fwo  the  day. 


SAGE  Volcanic  Aerosol 
Analysis  Methodology 
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Hampton  SAGE 
ABL  Range  Performance 
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Hampton  SAGE  range  performance  Is  similar  to  that  of  the  Middle  East.  This  implies  we 
can  use  the  LaRC  Lidar  data  to  fill  In  the  gaps  in  the  SAGE  data  for  both  locations. 
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LOS  ALAMOS  CONSULTING 
LOS  ALAMOS,  NM 


VIEWING  MODEL  USING  BLIRB  DATABASE 


WAVES  PROPAGATION  MODEL 
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23  km  Vis  (no  blurring)  10  km  Vis  (no  blurring) 


SPARTA’S  LIDAR  SIMULATION  CODE 
BACKSCAT  VERSION  4.0 
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nanweii  Avenue 
Lexington,  MA  02173 

Work  Performed  Under  Contract  FI  9628-C-9 1-0093 


WHAT  IS  BACKSCAT? 
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-  Accepts  Simulation  Inputs 

-  Saves  and  Recalls  Previous  Simulation  Conditions 

-  Displays  Results  Graphically 


BACKSCAT  APPLICATIONS 
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GROWTH  OF  BACKSCAT 


i  3 

I-I  ^  -£2 

r-  CO  T3 
C  ^  £ 

O  W  1  ^ 

CO  ^  O  D 

|i  9  I 

Is  3  I 

d  i 

O  LL  3 
U.  <  CQ 


4=  tr 

ro  W 

CO  0 
O  ^  ^ 

CM  ©  -D 

•i  I  s  s 

all? 


m  ro  o  < 
^  tP  m 


5  ^ 
o>  i 


o  m 

m  o 
?  8 


z  o  ^ 


■s 

TJ 

*5  w 
<  o 
«  w 

P 

O  © 
t3  < 

©  XJ 

©  g 

^  I 

u  O 

•8  i 

D  CO 
CO  3 


0) 

1  s 

O  (0 

O  3 

N 
•2  2 
u  (S 

»=  TJ 
UJ  3 

E  c 
B  2 

^  tr 


323 


-  SIgnal-to-Nolsa  Added 

-  Library  of  Detectors  Included 

-  Coherent  Doppler  Lidar  Simulated 

-  Water  Clouds  Included 

-  Estimates  of  Molecular  Absorption  Available 


NEW  FEATURES  IN 
BACKSCAT  4.0 
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NEW  FEATURES  IN 
BACKSCAT  4.0  (cont.) 
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Doppler  Analysis 


SNR  PERFORMANCE  MODEL 
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SNR  MODEL  ASSUMPTIONS 
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Local  Oscillator  Power  Is  Large  Enough  to  Provide 
Shot-Noise-Limited  Operation  of  the  Receiver* 

Unity  Mixing  Efficiency  of  Local  Oscillator  with  Return 
(*  Coherent  Doppler  Systems  Only) 


ESTIMATES  OF  RANGE  AND 
VELOCITY  ACCURACY 
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-  No  Frequency  Jitter  in  Local  Oscillator* 
(*  Coherent  Doppler  Lidar  Systems  Only) 


BACKSCAT  MAIN  MENU 
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Aerosol  Backscatter 
Raman  Scattering 
Coherent  Doppler 


BACKSCAT  OUTPUT  AND  PLOTS 
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Backscatter  vs.  Altitude 
Signal-to-Noise  vs.  Range 
Signal-to-Noise  vs.  Altitude 


FUTURE  UPGRADES 


POINTS  OF  CONTACT 
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(617)-377-3018 


THE  SOLAR  IRRADIANCE  BY  COMPUTATION 


Robert  L .  Kurucz 

Harvard  “Smithsonian  Center  for  Astrophysics 
60  Garden  St,  Cambridge ^  MA  02138 

I  am  now  able  to  compute  a  purely  theoretical  model  photo¬ 
sphere  (Kurucz  1992a;b;c)  that  reproduces  the  irradiance  measure¬ 
ments  of  Neckel  and  Labs  (1984)  in  the  visible  for  bandpasses  of 
approximately  2  nm.  That  model,  and  Avrett's  empirical  quiet 
sun  model  (Fontenla,  Avrett,  and  Loeser  1993)  that  includes  the 
chromosphere,  are  used  to  predict  the  irradiance  out  to  200 
microns  at  low  resolution. 

To  get  a  feel  for  the  scope  of  the  monochromatic  irradiance 
problem  I  have  computed  the  spectrum  from  150  nm  to  200'  microns 
at  a  resolution  of  500000  using  58  million  lines,  both  predicted 
and  observed.  If  this  spectrum  is  degraded  to  the  resolution  of 
the  model  it  looks  like  the  model.  At  any  given  wavelength  the 
spectrum  is  not  reliable.  But  at  a  resolution  of  10000,  say,  it 
approaches  measurement  accuracy.  In  regions  of  low  transmission 
it  is  more  reliable  than  existing  measurements.  I  will  publish 
tables  of  these  irradiance  spectra. 

I  am  producing  atlases  of  the  solar  flux,  central  intensity, 
and  limb  spectra  taken  by  James  Brault  at  Kitt  Peak.  One  atlas 
'•Solar  Flux  Atlas  from  294  to  13  00  nm"  by  Kurucz,  Furenlid,  Brault 
and  Testerman  (1984),  has  been  published  thus  far.  I  have  have  the 
Photometric  Atlas  of  the  Solar  Spectrum  from  1,850  to  10,100  cm-1** 
by  Delbouille,  Roland,  Brault,  and  Testerman  (1981)  also  taken  at 
Kitt  peak.  In  addition  I  have  the  ATMOS  central  intensity  atlas  from 
650  to  4800  cm-1  taken  by  Farmer  and  Norton  (1989)  from  Spacelab  3. 
The  replacement  for  the  flux  atlas  will  be  printed  on  demand  in  4 
parts  of  about  500  pages  each.  Each  page  will  show  the  observed 
spectrum  normalized  to  a  continuum,  the  state-of-the  art  computed 
transmitted  spectrum  on  the  day  the  atlas  is  printed,  and  line 
identifications.  Each  part  will  cost  $10  0,  plus  shipping  if  sent 
overseas.  If  I  can  obtain  a  color  laser  printer,  I  will  make  a 
$200  version  with  color  coded  solar  spectrum,  transmission  spectrum, 
transmitted  spectrum,  and  the  line  indent  if  icat  ions .  I  will  produce 
CD-ROMS  with  the  spectrum  and  line  data.  I  have  spectra  that  will 
continue  the  atlas  out  to  5  microns.  I  am  reducing  SMM  spectra  at 
shorter  wavelengths.  Similar  atlases  will  be  made  for  the  solar 
center  and  limb  in  collaboration  with  my  colleague  Barbara  Bell. 
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Parts  of  the  flux  atlases  directly  give  the  residual  irradiance 
spectrum  but  much  is  confused  or  obscured  by  terrestrial  lines.  I 
can  compute  them  away  except  where  the  transmission  is  too  low. 

Then  I  will  fill  in  with  the  purely  theoretical  solar  spectrum.  In 
this  way  I  will  finally  produce  an  atlas  showing  the  solar  spectrum 
above  the  atmosphere. 

I  use  these  atlases  to  test  the  pure  calculations  of  solar 
spectra  and  transmission  spectra.  I  identify  problems  with  the 
line  data  and  I  try  to  make  generic  corrections  that  improve 
hundreds  or  thousands  of  lines  at  a  time.  If  the  spectrum 
calculations  look  good  in  the  regions  of  high  transmission^  I 
can  have  some  confidence  that  the  regions  of  low  transmission 
are  computed  accurately.  The  main  problem  has  been  continuum 
placement.  Ozone  and  02  "dimer"  features  are  difficult  to 
determine  because  the  atlases  are  each  made  up  of  a  number  of 
sharply  peaked  FTS  scans.  The  continuum  placement  affects  the 
appearance  of  line  wings  and  the  apparent  depth  of  weak  features. 
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Ozone  cross-sections  for  the  Chappuis 
and  Wulf  absorption  bands 
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rovide  the  ozone  cross-  section  from  407  to  1089  nm. 

Temperature  dependence  is  also  given  for  407  to  762  nm. 


Spectral  Measurements  of  Ozone  Absorption 


Compilation  of  the  Data 
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Merging  overlapping  regions  with  smoothly  weighted  average 
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Highlights  of  Laboratory 
Measurements 


Ozone-friendly  single-pass  aluminum  absorption 
cell,  0.46m  opticd  path 

Fractional  absorption  detection  limits  below 
0.05%  (1  sec),  optical  depths  of  0-10% 

Visible  region 

Precision  capacitance  manometer  calibrated 
against  primary  pressure  standard  on-site 

Impurity  gas  pressures  measured  and  corrected 
for  in  re^-time 

Five  HeNe  laser  transitions  span  Chappuis  band, 
eliminate  wavelength  uncertainty 

Ozone  concentrations  cover  nearly  a  factor  of  20 
(  1.8  to  35  X  10*6  cm'^  as  well  as  zero) 

Near-IR  region 

Ozone  densities  determined  in  real-time  by 

absorption  at  632.8nm  HeNe  laser  transition 

Ozone  concentrations  cover  nearly  a  factor  of  10 
(  0.8  to  7  X  10^^  cm*^  as  well  as  zero) 
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Ratio  of  XS(T)/ 


T  dependence  of  03  vs  wavelength 


Temperature  [degK] 
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TABLE  1.  Absolute  Absorption  Cross  Sections  of  Carbon  Dioxide 
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Wave 


FAScode  for  the 
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Phillips  Lab/Optical  Physics  Division 
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Engineering 
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Spectral  Scanning  Using  FFT's 


New  Applications 
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Packaging 

•  All  Independent  Modules 
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August  31,  1993:  Pt.  Mugu,  CA 
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ALTITUDE  (KM) 


April  5,  1993:  State  College,  PA 


im)  3aniiiiv 


532  nm  BACKSCATTER  INTENSITY  607  nm  BACKSCATTER  INTENSITY 


ALTITUDE  (KM) 


I  t  1  I  I  I  I  t  I  t  I  I  i  -  1 . -  t  I  !  I  t  1 


CO 

CD 

LO 

a> 

O) 

05 

CJ) 

C75 

6 

d 

d 

d 

d 

(1AI9Z/)  NOISSIlAISNVdl 


368 


inatubo  Extinction  Profile  for  11/27/91 
4.8  degrees  north 
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Ground-Based  Measurements  of  HF  and  HCI 


Hilary  E.  Snell*  and  Paul  B.  Hays 

Space  Physics  Research  Laboratory 
Department  of  Atmospheric,  Oceanic,  and  Space  Sciences 
The  University  of  Michigan 
2455  Hayward 
Ann  Arbor,  Ml  48109-2143 

*Currently  with  AER,  Inc.,  Cambridge,  MA 


Abstract 

This  poster  describes  remote  sensing  observations  of  HF  and  HCI  column  densities 
over  Ann  Arbor,  Michigan.  Spectra  were  collected  between  October  1992  and 
July  1993  using  a  ground-based  Michelson  interferometer  and  the  technique  of  solar 
absorption  spectroscopy.  While  the  column  densities  computed  from  these  spectra 
agree  with  other  groups'  measurements  of  these  species,  we  have  noticed  that  the 
HCI  spectral  line  at  2925.897  cm'1  appears  strongly  asymmetric.  It  is  our  belief  that 
this  is  due  to  spectral  contamination  by  another  chemical  species  and  not  instrumental 
phase  errors  as  previously  reported  by  other  groups. 
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Introduction 


Hydrogen  chloride  (HCI)  and  hydrogen  fluoride  (HF)  are  important  gases  in 
stratospheric  chemistry.  Hydrogen  chloride  is  both  a  reservoir  and  a  sink  for  chlorine 
atoms,  which  play  an  integral  role  in  the  chemistry  of  ozone.  Hydrogen  fluoride  is 
chemically  inert  in  the  stratosphere  and  serves  as  a  sink  for  fluorine  atoms  through 
diffusion  to  the  tropopause  and  subsequent  rainout.  The  concentration  of  HF  is 
determined  by  the  destruction  of  anthropogenically-emitted  fluorine-containing 
compounds  while  the  HCI  concentration  has  both  anthropogenic  and  natural 
components.  The  relative  abundance  of  HF  and  HCI  provides  an  indication  of  the 
relative  importance  of  CFC's  to  other  sources  of  chlorine  in  the  total  amount  of  chlorine 
present  in  the  stratosphere. 

Remote  sensing  of  HCI  and  HF  has  been  routinely  accomplished  by  analysis  of 
the  near-infrared  spectrum  of  each  species'  fundamental  (1-0)  rotation-vibration  band. 
Numerous  lines  from  these  bands  have  been  used  for  atmospheric  density 
measurements  and  the  FASCODE  atmospheric  model  was  used  to  determine  the 
suitability  of  each  of  the  lines  for  ground-based  measurements.  Only  some  of  the  HF 
and  HCI  lines  can  be  used  in  ground-based  measurements  due  to  absorption  by 
methane  and  water  vapor.  It  is  for  this  reason  that  much  of  the  remote  sensing  of  HCI 
and  HF  has  been  done  from  balloons.  For  both  HF  and  HCI  the  Ri  transition  was 
selected  as  the  most  suitable  line  for  this  study,  2925.897  cm‘^  for  HCI  and 
4038.962  cm-1  for  HF. 
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Instrumentation 


The  measurements  described  in  this  poster  were  accomplished  using  a  Bomem 
DA-8  Michelson  interferometer  with  a  calcium  fluoride  (CaF2)  beamsplitter  and 
indium-antimonide  (InSb)  detector.  A  potassium-bromide  (KBr)  window  at  the 
instrument  entrance  aperture  enabled  instrument  evacuation  during  data  collection. 
An  optical  filter  was  placed  at  the  limiting  aperture  to  reduce  the  spectral  region 
incident  on  the  detector  and  allow  for  an  increase  in  the  amplifier  gain  settings.  Solar 
tracking  was  accomplished  with  a  roof-mounted  heliostat  system  and  the  light  was 
collimated  before  entering  the  instrument. 

The  interferograms  were  acquired  through  a  maximum  optical  path  difference  of 
200  cm  to  yield  an  unapodized  spectral  resolution  of  0.0025  cm’1.  This  optical  path 
difference  is  much  larger  than  the  interferogram-halfwidth  of  the  HF  and  HCI  features 
and  the  spectral  lines  are  completely  resolved.  Each  set  of  data  consists  of  two  co¬ 
added  interferograms  and  required  a  total  of  about  10  minutes  of  data  collection. 
Due  to  the  Bomem  software  and  computer  memory  limitations  only  data  from  a 
10  cm‘1  wide  spectral  region  was  saved. 

Spectra  representative  of  measurements  from  Ann  Arbor  are  shown  in  Figure  1 . 
The  spectrum  in  Figure  la  was  measured  on  1  October  1992  at  a  solar  zenith  angle  of 
67.04°  while  the  spectrum  in  Figure  1b  was  collected  on  20  January  1993  at  a  solar 
zenith  angle  of  78.26°.  As  indicated  the  HF  and  HCI  absorption  features  lie  on  the 
edge  of  very  strong  absorption  due  to  methane. 


Figure  1a  Spectrum  computed  from  the  interferogram  collected  at  8:44  est  on 
1  October  1 992.  This  is  the  result  of  two  scans  with  a  maximum  optical  path 
difference  of  200  cm,  which  yields  an  unapodized  spectral  resolution  of 
0.0025  cm-T 
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Figure  1b  Spectrum  computed  from  the  interferogram  collected  at  9:17  est  on 
20  January  1993.  This  is  the  result  of  two  scans  with  a  maximum  optical 
path  difference  of  200  cm,  which  yields  an  unapodized  spectral  resolution  of 
0.0025  cm-T 
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Determination  of  Column  Density 

The  column  density  was  retrieved  from  the  spectra  by  removing  the  background 
absorption  and  fitting  the  spectral  line  of  interest.  Removing  the  background  simplifies 
data  analysis  because  it  is  difficult  to  quantify  parameters  such  as  the  exact  field  of 
view  of  the  instrument,  detector  gain  and  spectral  response,  and  optical  losses  to  the 
signal  from  light-collecting  optics  and  within  the  instrument  itself. 

The  spectral  lines  were  fit  by  assuming  a  Voigt  lineshape  and  U.S.  Standard 
Atmosphere  temperature  and  pressure  profiles.  The  Air  Force  extension  to  the  U.S. 
Standard  Atmosphere  provides  the  HF  and  HCI  mixing  ratios  to  be  used  as  a  starting 
point  for  the  fit  to  the  spectral  data.  These  profiles  represent  a  mean  value  for  the 
desired  constituents.  The  fit  was  accomplished  by  allowing  variation  in  mixing  ratio 
magnitude  and  altitude.  Thus 

^measured[z]  =  M^af[z*(l  +  S)]  [1] 

where  ^af(z)  is  the  Air  Force  standard  profile,  M  accounts  for  a  magnitude  shift  in  this 
profile,  and  the  quantity  (1+5)  allows  for  vertical  transport  of  the  species  of  interest. 
Depending  on  the  value  of  5  the  mixing  ratio  profile  will  expand  or  contract  in  altitude; 
a  positive  value  of  5  indicates  subsidence  and  the  mixing  ratio  profile  will  be 
compressed.  This  method  of  varying  the  mixing  ratio  profile  is  more  realistic  than 
simply  shifting  the  profile  in  altitude.  Once  a  mixing  ratio  profile  has  been  determined 
it  is  used  to  compute  the  column  density. 
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Wavenumber  (cm'^) 

Figure  3  The  measured  HF  absorption  after  the  removal  of  the  background 
absorption  and  airmass  correction  factor  p. 
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Method  of  Removing  the  Background  Absorption 

The  measured  spectrum  may  be  written  as  the  product  of  background 
transmission  and  the  transmission  of  the  species  of  interest.  For  the  case  of  the 
spectral  regions  chosen  for  HF  and  HCI,  the  background  spectrum  consists  of 
absorption  due  to  methane  and  water  vapor.  Because  the  expression  for  the  Voigt 
lineshape  must  be  evaluated  numerically,  it  is  computationally  more  convenient  to 
implement  a  fitting  routine  using  the  analytic  solutions  for  the  Lorentz  and  Doppler 
shapes.  In  these  spectral  regions  the  water  and  methane  lines  are  strongly  absorbing 
and,  for  the  purpose  of  fitting  the  background,  we  have  treated  them  only  as  Lorentzian 
lines.  The  wings  of  these  features  will  still  have  a  Doppler  component  due  to 
stratospheric  absorption  but  exclusion  of  this  term  introduces  little  error  to  the  final  fit. 
The  initial  values  for  computation  of  the  Lorentz  shape  of  these  lines  are  tabulated  in 
the  HITRAN  database.  In  order  to  eliminate  the  HF  or  HCI  line  from  the  criteria  for  a 
good  fit  to  the  background,  only  the  spectrum  greater  than  three  halfwidths  away  from 
the  HCI  or  HF  absorption  was  included  in  the  background  fit. 

While  we  attempted  to  use  realistic  parameters  to  fit  the  background  we  soon 
found  that  doing  so  is  extremely  difficult.  In  particular,  we  found  that  (1)  the  methane 
linewidth  information  in  the  HITRAN  database  is  not  always  reliable;  (2)  line  mixing 
strongly  affects  the  wings  of  methane  lines,  particularly  in  the  region  around  the  HCI 
line  we  have  chosen;  and  (3)  line  mixing  is  most  severe  in  the  low-wavenumber  wings 
of  the  methane  P-branch  which,  again,  is  where  the  HCI  line  is  located.  To  overcome 
these  problems  we  decided  to  allow  the  half-widths,  magnitudes,  and  center 
wavenumbers  to  vary  while  computing  the  background  fit.  This  does  not  produce  a 
geophysically-meaningful  result  but  it  did  allow  us  to  remove  the  shape  of 
the  background. 
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Asymmetric  HCI  Line 

The  first  thing  one  notices  when  examining  Figure  5  is  that  the  HCI  line  is  not 
symmetric.  This  shape  has  been  noticed  by  other  groups  and  has  been  attributed  to 
instrumental  effects,  such  as  a  problem  with  the  phase  correction.  We  argue  that  if  the 
problem  were  instrumental  one  would  expect  the  other  lines  in  this  spectral  region  to 
be  distorted  as  well.  Unfortunately  this  spectral  region  does  not  contain  isolated  lines 
with  the  same  magnitude  of  absorption  as  the  HCI  feature.  Close  examination  of 
isolated,  but  stronger,  spectral  lines  indicates  that  the  other  features  in  this  spectrum 
are  not  asymmetric.  Examination  of  the  region  used  to  measure  HF  also  does  not 
show  any  obvious  asymmetries.  Furthermore,  this  shape  is  observed  in  all  of  our 
measurements  taken  over  several  months,  changes  magnitude  with  solar  zenith  angle 
in  the  manner  expected  for  a  spectral  feature,  and  appears  to  be  similar  to  spectra 
measured  by  other  groups  at  Mauna  Loa,  Hawaii,  and  Antarctica  (Rinsland  et  a!., 
1988;  Liu  ef  a/.,  1992). 

As  a  means  of  eliminating  instrumentation  errors  as  a  source  of  the  asymmetry, 
data  was  collected  using  different  mirror  scan  velocities  (0.5  cm  s""'  and  1  cm  s""*), 
different  beamsplitters  (CaF2  and  KBr),  different  widths  of  the  saved  spectral  region 
(2920-2930  cm‘1,  2922-2932  cm‘1,  and  2923-2929  cm‘1),  and  different  spectral 
resolutions  (0.0025  cm-1  and  0.005  cm""*).  In  all  cases  the  line  shape  is  asymmetric. 
Figure  6  illustrates  spectra  measured  with  the  KBr  and  CaF2  beamsplitters  at  the  two 
different  spectral  resolutions.  The  change  in  the  slope  of  the  lines  connecting  points 
clearly  indicates  that  the  spectral  line  is  asymmetric.  If  the  line  shape  was  symmetric 
these  connecting  lines  would  be  parallel;  the  lines  would  be  parallel  with  zero  slope  if 
one  of  the  resolution  elements  fell  exactly  on  the  peak  of  the  feature.  If  the  asymmetric 
shape  was  caused  by  a  problem  in  mirror  alignment  through  the  scan  one  might 
expect  the  asymmetry  to  be  reduced  if  the  measurement  is  made  at  a  lower  spectral 


resolution.  As  can  be  seen  in  Figure  6  this  is  not  the  case.  Thus  we  believe  that  the 
asymmetric  shape  of  the  line  is  real  and  is  not  an  artifact  of  the  instrument. 

in  order  to  support  the  hypothesis  of  the  presence  of  another  species,  the 
absorption  by  HCI  was  computed  and  removed  from  the  spectrum  shown  in  Figure  6. 
The  residual,  shown  in  Figure  7,  reveals  what  appears  to  be  another  spectral  line. 
Though  there  is  a  solar  CO  line  in  the  vicinity  of  the  HCI  absorption,  its  line  intensity  is 
several  orders  of  magnitude  weaker  than  that  of  HCI.  The  width  of  the  residual 
feature  is  approximately  0.047  cm’ "I  indicating  a  pressure  width  corresponding  to 
about  500  mb.  This  in  turn  indicates  that  the  gas  causing  the  absorption  is  uniformly 
distributed  with  altitude.  Furthermore,  because  it  is  present  in  spectra  obtained  at 
Antarctica  and  Mauna  Loa  it  is  not  a  chemical  species  endemic  to  Ann  Arbor. 
Unfortunately  we  have  been  unable  to  determine  exactly  what  species  accounts  for 
this  absorption  feature. 
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Figure  5  The  measured  HCI  absorption  after  removal  of  the  background 
absorption  and  airmass  correction  factor. 
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Figure  6  Comparison  of  data  collected  using  the  KBr  and  CaF2  beamsplitters  at 
unapodized  resolutions  of  0.0025  and  0.005  cm-^.  Background  absorption 
and  the  airmass  correction  factor  p  have  been  removed  from  the  spectra. 
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Figure  7  Difference  between  the  absorption  shown  in  Figure  5  and  the  computed 
HCI  absorption.  This  residual  illustrates  the  presence  of  another  spectral 
feature. 
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Retrieved  Densities 


utilizing  the  above  analysis  techniques  we  have  measured  the  column  content 
for  HF  and  HCI.  The  column  density  computed  from  each  of  the  spectra  collected  on  a 
given  day  was  averaged  and  the  (la)  standard  deviation  computed.  Information  about 
the  measurements  is  summarized  in  Table  1,  while  Table  2  summarizes  the  results  for 
HCI  and  Table  3  summarizes  the  HF  data. 

A  potential  source  of  error  in  determining  the  column  density  is  related  to  the 
calculation  of  the  solar  zenith  angle.  A  single  spectrum  is  derived  from  two 
interferograms  which  were  added  to  increase  the  SNR  and  requires  about  10  minutes 
to  record.  As  such  there  is  some  ambiguity  as  to  the  time  which  should  be  used  to 
compute  the  SZA.  To  estimate  the  magnitude  of  this  error  the  column  density 
computed  using  the  time  attached  to  the  data  file  was  compared  to  the  density 
computed  with  a  time  ±5  minutes  from  the  recorded  time.  Examination  of  this  error 
shows  that  for  most  cases  it  is  insignificant  compared  to  the  overall  scatter  in  the  data. 

The  number  of  spectra  collected  was  limited  by  several  factors.  As  one  might 
expect,  thick  clouds  make  it  difficult  to  collect  solar-absorption  spectra,  though  even 
near-sub-visual  cirrus  clouds  caused  problems  in  the  instrument  phase  correction 
routine.  Consequently  most  of  the  spectra  were  obtained  before  local  noon. 
Data  collection  was  further  limited  because  the  instrument  is  shared  by  several 
research  groups.  Furthermore,  measurement  of  HF  and  HCI  required  the  use  of 
different  optical  filters.  Changing  filters  requires  time  to  re-evacuate  the  instrument 
(a  sure  sign  that  it  would  soon  become  cloudy!).  Thus  on  any  given  day  data  was 
usually  obtained  only  for  a  single  spectral  region. 

Measurements  of  atmospheric  HCI  for  this  study  were  conducted  between 
October  1992  and  July  1993.  The  derived  column  densities  are  shown  in  Figure  10  as 
a  function  of  the  airmass  correction  factor.  Averaging  all  of  the  data  yields  a  column 
density  of  4.14±1.08  x  1015  molecules  cm'^.  The  large  standard  deviation  can  be 


attributed  to  errors  in  the  fit  due  to  the  asymmetric  shape  of  the  HCI  feature.  Errors  in 
tracking  due  to  the  variation  in  the  solar  zenith  angle  over  the  course  of  the 
measurement  were  found  to  be  negligible  compared  to  other  sources  of  error. 

Because  the  HF  line  is  symmetric  and  yields  a  much  better  fit  than  HCI, 
measurements  of  the  HF  column  density  were  made  to  a  higher  precision  than  the  HCI 
measurements.  The  data  from  Table  3  is  shown  in  Figure  11  as  a  function  of  the 
airmass  correction  factor.  Averaging  all  of  the  data  yields  a  column  density  of 
1.62±0.24  X  10^5  molecules  cm‘2.  The  standard  deviation  for  the  HF  data  is  much 
less  than  for  the  HCI  measurements  because  the  HF  line  is  symmetric  and  a  much 
better  fit  is  obtained.  As  with  the  HCI  measurements  the  errors  due  to  the  calculation 
of  the  solar  zenith  angle  are  insignificant. 


Date 


Day  Number 


HCI  Spectra 


HF  Spectra 


1 -OCT-1 992 

92-275 

43 

- 

27-DEC-1992 

92-362 

4 

7 

18-JAN-1993 

93-018 

- 

4 

19-JAN-1993 

93-019 

9 

20 

20-JAN-1993 

93-020 

- 

6 

1 -MAR-1 993 

93-060 

22 

- 

7-MAR-1993 

93-066 

- 

30 

18-MAR-1993 

93-077 

15 

- 

10-JUL-1993 

93-191 

26 

- 

13-JUL-1993 

93-194 

17 

Table  1  Summary  of  HF  and  HCI  spectra. 
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Day  Number 

Column  Density 
(molecules  cm- 
2) 

Estimated 

Accuracy 

Error  due  to 
SZA  Calculation 

92-275 

3.4412x1015 

19.42% 

1.75% 

92-362 

3.7538x1015 

25.58% 

0.38% 

93-019 

3.9644x1015 

15.81% 

3.55% 

93-060 

3.8136x1015 

11.94% 

3.65% 

93-077 

5.1441  x10i5 

9.05% 

3.52% 

93-191 

4.5841  x10i5 

36.23% 

2.51% 

93-194 

5.9122x1015 

1.90% 

2.81% 

Average  of  all 
data 

4.1395x1015 

26.06% 

2.74% 

Table  2  Summary  of  HCI  column  density  measurements. 


Day  Number 

Column  Density 
(molecules  cm- 
2)  _ 

Estimated 

Accuracy 

Error  due  to 
SZA  Calculation 

92-362 

1.2715  x  1015 

7.45% 

1.26% 

93-018 

1.3741  x  1015 

4.69% 

5.87% 

93-019 

1.5748  X  1015 

13.12% 

0.78% 

93-020 

1.6597x1015 

11.56% 

5.26% 

93-066 

1.7651  X  1015 

9.90% 

1.59% 

Average  of  ali 
data 

1.6212x1015 

15.10% 

1.84% 

Table  3  Summary  of  HF  column  density  measurements. 
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Comparison  with  Other  Measurements 

As  a  means  of  validating  the  ground-based  HF  and  HCI  measurements  made 
from  Ann  Arbor,  Michigan  we  have  examined  other  column  density  measurements  of 
HF  and  HCI.  These  comparisons  are  shown  in  Figure  12  for  HF  and  Figure  13  for  HCI. 

One  source  of  recent  data  for  comparison  is  the  Halogen  Occultation 
Experiment  (HALOE).  One  of  several  instruments  aboard  the  Upper  Atmosphere 
Research  Satellite  (UARS),  HALOE  uses  solar  occultation  techniques  to  measure  the 
atmospheric  mixing  ratio  profiles  of  several  gas  species,  including  HF  and  HCI.  The 
instrumental  requirements,  experimental  objectives,  and  geographical  coverage  are 
described  by  Russell  et  al.  (1993).  We  have  examined  HALOE  mixing  ratio  profiles 
measured  close  to  the  location  of  Ann  Arbor  (42.28°  N,  83.71°  W).  The  lower  altitude 
limit  for  the  HALOE  data  is  about  14  km  and  the  mixing  ratio  profile  can  be  integrated 
to  yield  a  column  density  from  this  point  to  the  top  of  the  atmosphere.  From  this  data 
the  column  density  for  HF  is  1.76x10"*  ^  cm‘2  while  the  HCI  column  density  is 
2.23x10^5  cm'2.  These  values  are  accurate  to  10-15%  (Russell  et  al.,  1993). 
Hydrogen  fluoride  does  not  have  a  significant  tropospheric  source  and  the  HALOE 
value  falls  well  within  the  error  bars  of  our  measurement  (1.62±0.24x  10"*^  cm'2). 
Hydrogen  chloride  does  have  a  tropospheric  component.  If  the  density  decreases  at 
the  rate  projected  by  the  AFGL  atmospheric  profile  the  HCI  column  density  we  have 
measured  (4.14±1.08  x  10"*^  cm"2)  is  actually  too  low  when  compared  to  HALOE 
data  (see  the  Air  Force  column  density  profile  shape  in  Figure  13). 

The  other  comparisons  shown  in  Figures  12  and  13  illustrate  that  the  column 
densities  we  have  measured  are  comparable  to  previous  measurements.  However, 
due  to  differences  in  latitude,  season,  and  year  we  would  not  expect  these  values  to 
be  in  exact  agreement. 
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Because  the  HF  line  is  symmetric  and  yields  a  much  better  fit  than  HCI,  measurements 
of  the  HF  column  density  were  made  to  a  higher  precision  than  the  HCI  measurements.  The 
data  is  shown  as  a  function  of  the  airmass  correction  factor.  Averaging  all  of  the  data  yields  a 
column  density  of  1 .62±0.24  x  10i5  molecules  cm-2.  The  standard  deviation  for  the  HF  data  is 
much  less  than  for  the  HCI  measurements  because  the  HF  line  is  symmetric  and  a  much 
better  fit  is  obtained.  As  with  the  HCI  measurements  the  errors  due  to  the  calculation  of  the 
solar  zenith  angle  are  insignificant. 
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Column  Density  (molecules  cm  "2) 


Measurements  of  atmospheric  HCI  for  this  study  were  conducted  between  October 
1992  and  July  1993.  The  derived  column  densities  are  shown  as  a  function  of  the  airmass 
correction  factor.  Averaging  all  of  the  data  yields  a  column  density  of 
4.14+1.08  X  lO^s  molecules  cm-2.  The  large  standard  deviation  can  be  attributed  to  errors 
in  the  fit  due  to  the  asymmetric  shape  of  the  HCI  feature.  Errors  in  tracking  due  to  the 
variation  in  the  solar  zenith  angle  over  the  course  of  the  measurement  were  found  to  be 
negligible  compared  to  other  sources  of  error. 
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Summary 

Using  the  technique  of  solar  absorption  spectroscopy  we  have  measured  the 
atmospheric  absorption  due  to  HF  and  HCI.  From  this  information  we  are  able  to 
calculate  the  column  density.  Comparison  of  our  data  with  that  of  other  groups 
indicates  that  the  values  are  in  good  agreement  with  other  measurements  at 
this  latitude. 

Examination  of  the  HCI  absorption  feature  at  2925.897  cm'"'  under  various 
experimental  conditions  indicates  that  the  observed  asymmetry  in  the  shape  of  the  line 
is  not  introduced  by  the  interferometer.  Determination  of  the  species  causing  this 
absorption  is  necessary  if  this  spectral  line  is  to  be  used  to  measure  the  density  of  HCI. 
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AN  ATMOSPHERIC  MODEL  FOR 
GRAVITY  WAVE  INDUCED 
TURBULENT  LAYERS  (BLINI) 
BASED  ON  THE 

SATURATED  CASCADE  MODEL 


E.  DEW  AN  PLyGPOS 

N.  GROSSBARD  BOSTON  COLLEGE 
T.  VANZANDT  NOAA 


EDE6134GPO 


394 


OUTLINE 


•  PROBLEM  TO  BE  SOLVED 

•  SOLUTION  TO  PROBLEM 

•  GRAVITY  WAVE  SIMULATIONS 

•  RESULTS 

•  CONCLUSIONS 
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PROBLEM  TO  BE  SOLVED 


•  AF  SYSTEMS  USING  HIGH  POWERED  LASER 
BEAMS  ARE  LIMITED  BY  ATMOSPHERIC 
TURBULENCE 

•  TURBULENCE  OCCURS  IN  THIN  LAYERS  CALLED 
"BLINI"  (RUSSIAN  FOR  "PANCAKES”) 

DIMENSIONS:  100  M  VERTICAL 

100  KM  HORIZONTAL 

•  CAUSE  OF  TURBULENCE  LAYERS  =  GRAVITY 
WAVES  [BY  HYPOTHESIS] 

•  PROBLEM:  HOWTO  SIMULATE  "BLINI"  BEHAVIOR? 
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SOLUTION  TO  PROBLEM 


•  LINEAR  INSTABILITY  DUE  TO  WAVE  SHEARS 

Ri.^^0.5 

S 

•  USE  SATURATED  CASCADE  THEORY  TO  SIMULATE 
WAVE  FIELD 

•  REALIZATIONS  MUST  BE  BASED  ON  THE 
THEORETICAL  SPECTRA 

•  THE  SHEAR  REGION  WHERE  S'  >  2N'  WILL  DEFINE 
THE  BLINI  LOCATION  IN  X-Z  (DISTANCE-ALTITUDE) 
SPACE  [N  =  CONST.] 
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SATURATED  CASCADE  SPECTRA  AND 


'Pv,  (kz)  =  a  kz=‘  ■  27c  =  (5  X  1 0“*)  k? 

Vx  =  HORIZONTAL  VEL.  FLUCT. 

1^^  1^3 1^2^^  ^2501^ 

Va,N  J _ 

e  =  2x10'' mV"  (STRATOSPHERE),  |f=15.9 
N  =  2  X  1 0'"  s’  (STRATOSPHERE), 
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REALIZATIONS  OF  v,  FROM  (kj 

- -  -  ^ 


M 

v:  =  I  V,  (k.  (I)) 

1  =  1 

Vx  (kz)  =  A  (kz)  •  sin  (kx  X  +  kz  Z) 


(ALT.) 


(NOT  FOR  ENDPOINTS) 
L  =  TOTAL  LENGTH  =  N'AZ  (AZ  =  SPACING) 


l[kjl)  +  k,(l  +  i)] 


A'(k.(l))-J  K(K)dk,] 


©[k.(i)  +  k.(i-i)] 
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BLINI  OCCUR  IF: 


dv; 

dz 


(x,2)  >  2.8  X  1 0-'  s-' 


ALTERNATIVE  APPROACH 


M 


S  (X,Z)  =  2  B(k^  (I))  sin  (k^X  +k^Z) 

i  =  1 


f[k,(i)  +  k^(l  +  1)] 
(|)[kz(l)  +  k,(l-i)] 


TO  PROVIDE  RANDOM  PHASES  AND  AMPLITUDES. 
MULTIPLICATION  OF  A  (k,)  OR  B  (kj  BY  FOURIER 
TRANSFORMED  WHITE  GAUSSIAN  NOISE  IS  USED. 
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THE  FOLLOWING  SIMULATIONS  USE  THE 
SATURATED-CASCADE  SPECTRAL  MODEL 

a-  X-Y  PICTURE  OF  BLINI 
b-  SHEAR  SPECTRUM 


VERTICAL  SPACING  =  -^  =  600  m  ,  X,  =  100  m 

HORIZONTAL  SPACING  =  ^  =  285  km  ,  \  =  100  m 
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130  k’s  RANDOM  PHASES  AND  AMPLITUDES 


wavEN'JMEER  ( : .  /ri  ] 
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X  (KILOMETERS) 

200  k's  RANDOM  PHASE  AND  AMPLITUDE 
(VELOCITY) 
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IN'TEG^^PL  or  PSCr  1.24E-04 
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CONCLUSIONS 


•  WE  HAVE  EXPLORED  THE  BEHAVIOR  OF  TURBULENT 
BLINI  THAT  ARE  DUE  TO  GRAVITY  WAVES. 

•  THE  SATURATED  CASCADE  THEORY  WAS  USED  TO 
PRODUCE  SIMULATIONS. 

•  BLINI  APPEAR  TO  BE  PERIODIC  AND  BLINK  ON  AND  OFF 
WITH  RESPECT  TO  X  (OR  t)  AS  WELL  AS  WITH  ALTITUDE: 

SPACING  IS  AT  V2  and  ^2  DISTANCES 
RESPECTIVELY. 

•  THEY  SLOPED  FROM  UPPER  LEFT  TO  LOWER  RIGHT 
WITH  OUR  MODEL. 

•  THESE  SLOPES  WERE  CONSISTENT  WITH  THE  MOST 
DOMINANT  FREQUENCIES  IN  THE  SHEAR  SPECTRA. 

•  IT  WAS  FOUND  THAT  BLINI  OCCURRED  WITH  AND 
WITHOUT  THE  USE  OF  RANDOM  AMPLITUDES  IN  THE 
PSD'S  USED  TO  PRODUCE  THE  REALIZATIONS.  THIS 
WAS  NOT  EXPECTED.  IT  DID  NOT  OCCUR  IN  THIS 
MANNER  WITH  A  LESS  SOPHISTICATED  MODEL  USED 
EARLIER  (THE  “SEPARABLE  MODEL"). 

•  MULTIPLE  PERIODS  WERE  FOUND  IN  BLINI  WITH  LOW 
NUMBER  OF  WAVE  NUMBER  SIMULATIONS  AND  THEY 
WERE  CONSISTENT  WITH  DOMINANT  FREQUENCIES  IN 
THE  SHEAR  SPECTRA. 

•  EVIDENCE  FOR  THE  SATURATED-CASCADE  THEORY. 
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ALTITUDE  (  km  ) 


giai 


SlijSIS 


iUx^X^ 


s'^41 


i»^  •■•  ■'  ‘■-..‘•il 


LOCAL  TIME  (AST) 

WOO'D MnN  ,'RASTOGI_,  /?xii> Sato (W 

Figure  1.  Radar  echo-power  of  backscatter  signals  from 
turbulent  fluctuations  of  clear  air  in  the  stratosphere 
and  upper  troposphere.  Shade  levels  are  every  4  dB.  A 
piecewise  linear  trend  has  been  subtracted  with  0,  2,  5, 
and  12  dB  of  attenuation  at  9,  14,  19,  and  24  km,  respec 
tively.  Results  obtained  with  the  430  NHz  Arecibo  radar 
at ^iSHjneter  resolution.  Unpublished  material  courtesy 
of  Woodman,  Rastogi  and  Sato. 
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ARC  is  1-D  model,  estimate  effects  for  various  realistic  cases  (scale  sizes  of 
cloud/clear  regions 

Applicability  to  other  emitters 


Fig.  3.  Integrated  radiance  from  the  1-0  fundamental  bands 
or^-C'^0.  and  as  determined  from  spectral 

fits  to  numerous  scans  from  the  CIRRIS  lA  database,  under 
nighttime  (top)  and  daytime  (bottom)  conditions.  The  lines 
superimposed  on  the  data  points  are  predictions  made  by  the 
Atmospheric  Radiance  Code  (ARC)  model  [Wimersteiner  et 
aL.  1992]. 
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Limb  Line-ol-Sight  Radiance 


Normalized  Radiance  Radiant  Enhancement 
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-  CO  mixing  ratio  profiles  must  be  supplied  to  model 

-All  CO  profiles  exhibit  unusual  behavior,  with  large  Increase  in  mixing 
ratio  in  stratosphere  and  mesosphere 

-  Large  uncertainty  in  CO  mixing  ratio  profiles 
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NIGHT  / 


VIBR.\TIONAL  TEMPERATURE  (K) 

Fig,  4.  Nighttime  (top)  and  daytime  (bottom)  vibrational 
temperature  profiles  for  and  as 

predicted  by  the  ARC  model.  Also  shown  for  reference  is 
the  kinetic  temperature  input  to  the  ARC  calculation. 
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C’®0  Sensitivity  to  Lower  Boundary 


(uj)i)  epnjijiv 
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Temperature  (K) 


t3C^®0  Dependence  on  Cloud  Structure 
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Temperature(K) 


Limb  Band  Radiance  of  '^C '°0  and  '®0:  July  10 
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edge  effects 

-  Structured  radiance  from  lower  boundary  clouds  likely  to  be  at  least  as 
important  (nighttime)  as  in  situ  variability  (gravity  waves  ->  T,  [M]) 


INFRARED  RADIANCE  FLUCTUATIONS 
IN  THE  UPPER  ATMOSPHERE 


I — 

CC  LU 
LU  Q 
CD  _ 

DC  < 
-  Q_ 

£  Q 
CD  Z 
Z  < 


LU  O 


>-  ^ 
DC 

O  cn 

h-  LL 

<  < 

CC  ^ 

o  ^ 

DQ  O 
<  CJ 
-J  CD 

O-  < 
“i  I 


448 


JUNE  8,  1994 


(f) 

LU 

to 

CO 

< 

CL 

Q 


O 

< 

o 

DC 

CL 

CL 

< 

LL 

O 


in 

y 

I — 

y 

h- 

< 

h- 

in 

Q 

< 

in 

z: 

o 

H 

< 

h" 

DC  O 
H  3 


LU 

o 


^  I 

LU  ^ 

-  LU 

CO  >  X 

^  CE  ^  < 

<  UJ  ^  5 

O  >  < 

CO  O  ^ 


< 

DQ 

yi 

DC 

DC 

O 

CO 


< 

h- 

< 

Q 


LU 

Q 

O 


LU 

Q 

O 


O  O  ^ 

I — 

< 


h- 

X 

y 

LU 


in 

in 


O  $ 

_ I  LL 

<  o 
O  ^ 

LU  O 
CL  O 

5::  Q 

LU 
DC 
0- 


X 

H- 


in 


in 

LU 

o 

< 


in 

in 

< 

0. 

Q 


O 

y 

DC 

<  _ 

CL  □_ 


< 

GD 

LU 


>- 

DC 

< 


< 

CO 


o  < 

u  in 


CO 


449 


C/D 

U 

h- 

(f) 


C/D 

LU 


H 

< 

h- 


LU 


Q 
D 
I — 
C/D 


Q 
Z 
< 

^  m 

o  >- 
O  CO 

I-  > 

H  <  c/5 

y  25 

gom 

DC/D 
rrJ  D 


DC 

CL 


O  LU 

H-  U 


> 
I — 


DC 

o 

LU 


DQ 

< 

CL 

< 

O  CC 

LU 


CO 
LU 

Q  CD 


< 

d: 


o  — 

DC 


LU 

X 

H- 

CL 

O 

_I 

LU 


X 

CL 

C/D 

O 


LU 

h- 

< 

cc 

LU 


h-  UJ 

>  <  o 

LU  LU  O 

Q  O  H 


CO 


z  s 

<  o 

?  3 

U  LU 

3  i 
t  5 


< 

DC 

LU 

D- 

LU 

I- 

co 

< 

CD 

D 

CL 


H 

< 

DC 

DQ 

> 

I 

o 

DC 


Q 

CO 

Q_ 

Q 


CO 

Z 

LU 

CD 
u 
< 

I- 

O 

w 

LU  _J 
DC  < 

b  o 


D 
O. 

CL  < 
LU  LU 

H  O 
< 

H 


< 

DC 

< 

> 

o 

CJ 


LU 

O 

z 

< 

Q 

< 

DC 


H  I — 
O  O 

Q  Q 

LU  LU 
DC  DC 
CL  0. 


LU 

H 
< 
1 — 
CO 


>- 

Q 


< 

LL! 

H 

LU 

CO 

h- 

CQ 

X 

h- 

o 

z 

LU 

_J 

LU 

o 

DC 

LU 

H 

o 

DC 

Q 

LU 

i- 

D 

CO 

LU 

u 

LU 

1 — 

O 

DC 

< 

o 

CO 

H 

< 

LU 

1“ 

LU 

A 

A 

CO 

LU 

LU 

LU 

DC 

h- 

_J 

— J 

H 

LU 

< 

< 

o 

DC 

O 

O 

z 

D 

CO 

CO 

CO 

o 

1 — 

LU 

LU 

U 

u 

1 

CC 

1- 

1 

i- 

1 

H" 

I 

CO 

r 


450 


CONSTRUCT  IMAGE  REALIZATIONS  OF  STRUCTURE 
-ARBITRARY  BAND  PASS  AND  LOS 


SHARC  APPROACHES 
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TODAY'S  TOPIC  IS  LOS  WEIGHTING 


OVERVIEW  OF  APPROACH 
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VIBRATIONAL  TEMPERATURE  DERIVATIVES 
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ALTITUDE  (KM)  ALTITUDE  (KM) 


OVERVIEW  LOS  RADIANCE  STATISTICS 
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SCALING  OF  "OVERLAYS"  IS  INADEQUATE 


FIELD  OF  VIEW  EFFECTS  ON 
CORRELATION  LENGTHS 


)RRELAT 
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INITIAL  MODEL  APPLICATIONS 
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ANDS 


CIRRIS-1A  RADIOMETER  BANDPASSES 


25%  spectral  response  limits 
USED  IN  CURRENT  SIMULATIONS 


HORIZONTAL  RADIANCE  CORRELATION  LENGTHS 

CIRRIS-1A  (4.1  -  4.5  //m) 
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CORRELATION  LENGTH  (KM) 


LOS  WEIGHT 
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VERTICAL  RADIANCE  CORRELATION  LENGTHS 

CIRRIS-1A  (4.1  -  4.5 //m) 


460 


CORRELATION  LENGTH  (KM) 


RADIANCE  STANDARD  DEVIATION 

CIRRIS-1A  (4.1  -  4.5 //m) 


PERCENT  (Jl/L 


HORIZONTAL  DISTANCE  (km)  §12 


HORIZONTAL  RADIANCE  CORRELATION  LENGTHS 

CIRRIS-1A  (2.6  -  3.3 //m) 
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CORRELATION  LENGTH  (KM) 


HORIZONTAL  RADIANCE  CORRELATION  LENGTHS 

CIRRIS-1 A  (8.0  -11.8  //m) 
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CORRELATION  LENGTH  (KM) 


HORIZONTAL  RADIANCE  CORRELATION  LENGTHS 

CIRRIS-1A  (8.5  -  18//m) 
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CORRELATION  LENGTH  (KM) 
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-  MUCH  DATA  WITH  NO  PREDICTIONS 


CONCLUSION 


in 

O 


in 

in  Q- 

P  ? 

<  CO 

H  < 

^  if) 
LU  O 
O  H 
Z  ^ 

% 

or  LU 

in  3 

z  < 

^  cc 

^  LU 

CC  % 

LU  ^ 

h- 

LU  r— 

Q  C/) 
LU 


I- 

co  2 

LU  LU 

z  Q 

LU  Z 

o  m 

CO  ^ 

DC  Q 
(T  CO 

o  2 

It  ^ 

3  i 

H  < 
D  “ 
O  Q 


Q  > 
O 


to  O 

9  I 

>  z 

O  LU 

d:  o 

0.  CO 


o  I 

LU  “ 

i  ° 

^  Q 
DC  LU 
LU  CQ 

s? 


1 1 1  LL 

O  LU 

CO  _ 
CO  LU  X 

I  W  ^ 

b  S  S 

^  CO  LJ 


>- 

< 

> 

X 

< 

O 

CO 

CO 

< 

u 

1- 

X 

z 

> 

o 

o 

CO 

LU 

CO 

iD 

LU 

X 

LU 

Q 

3 

H 

z 

o 

K 

CO 

< 

X 

Cl 

1— 

z 

Z 

z 

i- 

_J 

LU 

Q 

< 

CQ 

O 

< 

Z 

LU 

o 

X 

H- 

CL. 

H 

H 

Z 

LU 

LU 

Z 

< 

lj 

> 

LU 

z 

> 

1- 

CO 

z 

—1 

X 

CD 

co 

z 

< 

o 

LU 

1 — 

Q 

X 

CO 

1  I 

I  t 


CO  ^ 

^  z 
O  < 

b  to 

<  CO 

O  > 

□  ^ 

^  < 


467 
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